OFFICE 

OF 

RESEARCH 


iiW;  llli 


LEHIGH  UNIVERSITY  \M|) 


tv*  nmn  g’rs 

w 


THE  FATIGUE  STRENGTH 
OF  FILLET  WELDED  CONNECTIONS 


Karl  H.  Frank 


October  1971 


Fritz  Engineering  Laboratory  Report  No.  358.37 


NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

St  >  r.  ;fi«'  d  V.i  ,’2|M 


_ Unclassified _ 

SciunU  I'hissituMtum 

DOCUMENT  CONTROL  DATA  ■  R  &  D  '  ” 

*  .irifi  «  /.i'*  -.i/ii  •  * f  tttlw  luul\  nf  ,th’.;r,n  t  «*/!</  iiiilruiv  .iimuMfi  >  >n  nm  -I  e/ifrrit/  u/irn  thr  iMrMlJ  rr/j.»fr  j  .  t  /,is.w/h’«/| 


Lehigh  University 


*'*l.  Kl  *  l»K  t  SI  CUldh  (  l  A  •!  I  H  A  I  I 

Unclassified 


h  *  M  O  »  «  i  *  l  I 


THE  FATIGUE  STRENGTH  OF  FILLET  WELDED  CONNECTIONS 


4  I  s *.  n  i »  list  NO  M  s  |  T»  /»«•  «*f  r i*|*i *rf  ■mil  m«  fii  1 1 «  «/.< f «•  s  » 


Ak»  r"OWlsi  i  Firs  f  riiifpc,  nintdl*'  mi  fr«i/.  f/i  s  f  ri«irnt> 


Karl  H.  Frank 


Nt^ONT  ;  All 


»  O  '  A  l  '4i  ’  f  i  »  I  A  i,  t  \. 


».  )  O  *  III  I  ' 


-.1  ,N  n,*(  t  OK  bNAS’  NO 

N  000 14-68- A- 0514 

h  r  W  OjI  C  r  NO 


.  »  ‘  !•  '  VM  l.  ,’l 


Fritz  Eng.  Lab.  Rep.  No.  358.37 


if  Dim  »<  >«i  con  '  ’iOi'jI  i4/ii  nf/it-r  mifii/HTs  fh.if  n;.n  /••■  ti  *.  t&ini 
thi'-  rr/i.-rfj 


S  T  KIM  '  j  TION  STA’rigif^T 


"Approved  for  public  release;  Distribution  unlimited". 


1  *  I  *  A  M  •  *.  *  | 


Office  of  Naval  Research 


This  report  describes  a  study  of  the  fatigue  behavior  of  fillet 
welded  joints  stressed  perpendicular  to  the  weld  line.  The  study 
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behavior  of  load  and  non-load  carrying  fillet  wold  joints  was  determined 
This  experimental  study  concentrated  on  the  fatigue  behavior  in  the  r 
transition  region  between  high  and  low  cycle  fatigue  (10^  cycles  to  1 0 3 
cycles  of  load  application) . 
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ABSTRACT 


This  report  describes  a  study  of  the  fatigue 
behavior  of  fillet  welded  joints  stressed  perpendicular  to 
the  weld  line.  The  study  included  an  experimental  phase  in 
which  the  stiess-life  and  cracking  behavior  of  load  and  non¬ 
load  carrying  fillet  weld  joints  was  determined.  This 
experimental  study  concentrated  on  the  fatigue  behavior  in 

the  transition  region  between  high  and  low  cycle  fatigue 
3  5 

(10  cycles  to  10  cycles  of  load  application). 

The  second  phase  of  the  study  was  the  determination 
of  the  fracture  mechanics  stress  intensity  factor  for  cracks 
at  the  weld  toe  or  root  of  fillet  welded  joints.  The  finite 
element  technique  was  used  to  determine  the  compliance  of 
typical  fillet  welded  joints.  The  results  of  this  analysis 
showed  that  the  stress  intensity  factor  was  a  function  of 
weld  size,  weld  angle,  plate  thickness  and  weld  penetration. 
Relationships  to  estimate  the  stress  intensity  factor  were 
developed . 


The  third  phase  of  the  study  was  the  correlation 
of  the  experimental  fatigue  stress-life  behavior  of  fillet 
welded  joints  with  the  predicted  behavior  using  the  concepts 


r  ' 


-1 


ABSTRACT 


This  report  describes  a  study  of  the  fatigue 
behavior  of  fillet  welded  joints  stressed  perpendicular  to 
the  weld  line.  The  study  included  an  experimental  phase  in 
which  the  stress-life  and  cracking  behavior  of  load  and  non¬ 
load  carrying  fillet  weld  joints  was  determined.  This 
experimental  study  concentrated  on  the  fatigue  behavior  in 

the  transition  region  between  high  and  low  cycle  fatigue 
3  5 

(10  cycles  to  10'  cycl  ■'s  of  load  application)  . 

The  second  phase  of  the  study  was  the  determination 
of  the  fracture  mechanics  stress  intensity  factor  for  cracks 
at  the  weld  toe  or  root  of  fillet  welded  joints.  The  finite 
element  technique  was  used  to  determine  the  compliance  of 
typical  fillet  welded  joints.  The  results  of  this  analysis 
showed  that  the  stress  intensity  factor  was  a  function  of 
weld  size,  weld  angle,  plate  thickness  and  weld  penetration. 
Relationships  to  estimate  the  stress  intensity  factor  were 
developed . 


The  third  phase  of  the  study  was  the  correlation 
of  the  experimental  fatigue  stress-life  behavior  of  fillet 
welded  joints  with  t-he  predicted  behavior  using  the  concepts 
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of  the  fracture  mechanics  of  stable  crack  growth.  The 

study  showed  that  the  fatigue  behavior  of  various  joint 

geometries  could  be  correlated  using  the  stress  intensity 

factors  developed  in  phase  two.  The  correlation  included 

failures  from  cracking  at  the  root  and  toe  of  the  fillet 

weld.  The  analysis  provided  a  good  estimate  of  the  observed 

3  7 

fatigue  behavior  for  fatigue  lives  from  10  to  2x10  cycles 
of  load  application. 
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1 .  INTRODUCTION 


This  dissert ition  is  concerned  with  the  study  of 
the  fatigue  behavior  of  fillet  welded  joints.  Historically 
this  subject  has  held  the  interest  of  many  investigators. 
Interest  was  generated  by  the  general  use  of  fillet  welded 
joints  in  welded  construction.  Fillet  welds  are  usually 
more  economical  than  other  types  of  welds  and  are  exten¬ 
sively  used.  However,  experimental  studies  showed  that 
fillet  welded  joints  had  low  fatigue  strengths  when  they 
were  stressed  in  a  direction  perpendicular  to  the  weld 
line.  This  led  to  further  studies  in  attempts  to  determine 
why  these  reductions  were  observed. 

Early  studios  of  the  fatigue  behavior  of  fillet 

welded  joints  showed  that  their  fatigue  lives  were  a 

(1-7) 

function  of  weld  size  and  joint  geometry.  Photoelastic 

studies  were  used  to  evaluate  the  influence  of  the  joint 

(8  9 ) 

geometry  upon  the  state  of  stress  in  the  joint.  '  These 
studies  provided  insight  into  the  stress  distribution 
within  the  joints  and  provided  an  indication  of  the  stress 
concentrations  which  exist  at  the  toe  and  root  of  the  weld. 
Nearly  all  of  these  stress  studies  were  concerned  with  the 
static  behavior  of  the  joint.  They  did  not  provide  a 
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relationship  between  joint  geometry  and  applied  stress  and 
their  relationship  to  fatigue  life. 

The  increased  use  of  welded  construction  during 

the  late  1950's  focused  attention  on  the  fatigue  behavior 

of  fillet  welded  joints.  Numerous  investigations  were 

undertaken  on  different  joint  geometries,  steels  and 

welding  processes.  A  review  of  these  studies  is 

given  in  Kef.  26.  This  work  showed  that  the  principle 

variables  influencing  the  fatigue  life  of  a  fillet  welded 

joint  were  the  geometry  of  the  joint  and  the  applied  stress 

range.  The  static  strength  of  the  base  metal  was  observed 

to  have  little  influence  on  the  fatigue  strength  for  the 

range  of  steels  used  in  bridge  and  building  construe- 
(2526) 

tion.  '  '  ;  In  addition,  the  type  of  welding  process 

employed  in  the  manufacturing  of  the  joint  as  well  as  the 
strength  of  the  deposited  weld  had  little  influence.  They 
wei  observed  to  only  influence  the  fatigue  behavior  if  they 
chaf  jd  the  weld  profile  significantly  and/or  provided  a 
weld  with  different  sized  defects  . 

Recently  fracture  mechanics  has  been  employed  to 
analyze  the  fatigue  behavior  of  fillet  welded 
joints. (25,28,29 , 30)  These  studies  yielded  qualitative 
information  about  the  influence  of  joint  geometry  upon  the 
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observed  fatigue  behavior.  However,  adequate  estimates  of 
the  fracture  mechanics  stress  intensity  factor  for  cracks 
in  fillet  welded  joints  were  not  available.  Quantitative 
estimates  of  fatigue  life  behavior  require  more  accurate 
estimates  of  the  stress  intensity  relationship. 

The  results  of  these  previous  studies  have  not 
provided  a  means  to  quantitatively  estimate  the  fatigue  life 
of  fillet  welded  joints.  Each  set  of  test  results  must  be 
treated  separately  since  no  means  exists  to  correlate  the 
behavior  of  joints  with  different  geometries.  The  influence 
of  weld  size  and  penetration,  plate  thickness  and  initial 
flaw  size  has  not  been  studied  analytically.  This  disser¬ 
tation  describes  the  results  of  an  analytic  and  experimen¬ 
tal  study  of  the  fatigue  behavior  of  two  simple  fillet 
welded  joints  stressed  perpendicularly  to  the  weld  line.  A 
method  of  quantitatively  estimating  the  influence  of  joint 
geometry  upon  the  fatigue  life  behavior  of  fillet  welded 
cruciform  joints  is  developed.  The  study  consisted  of 
three  major  phases. 

The  first  phase  was  an  experimental  investigation 
to  determine  the  low  cycle  fatigue  behavior  of  two  fillet 

welded  joints.  Previous  studies  examined  the  high  cycle 
life- low  stress  region.  This  study  examined  the  transition 
region  between  low  cycle  life  (less  than  1000  cycles)  ana 
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high  cycle  life  (greater  than  100,000  cycles)  of  high 
strength  steel  fillet  welded  joints.  The  purpose  of  this 
phase  of  the  study  was  to  determine  the  stress-life 
behavior  and  the  mode  of  failure. 

The  second  phase  of  the  study  involved  a  stress 
analysis  of  the  test  joints  using  the  finite  element 
technique.  The  purpose  was  to  examine  the  effect  of  joint 
geometry  upon  the  stress  distribution  in  the  joint  and 
provide  fracture  mechanics  stress  intensity  factors  for  the 
cracking  behavior. 

The  third  phase  of  this  study  correlated  the 
fatigue  behavior  of  the  test  joints  with  their  geometry. 
This  was  accomplished  by  relating  the  stress  intensity  fac¬ 
tors  from  the  finite  element  study  to  the  crack  growth  rate. 
The  relationships  between  fatigue  behavior  and  joint 
geometry  permitted  the  results  of  other  studies  to  be 
correlated.  In  the  case  of  certain  welded  joints,  the 
influence  of  geometry  upon  the  fatigue  can  be  estimated 
quantitatively.  Clearly  weld  imperfections,  such  as  weld 
border  cracks  and  weak  regions  along  solidification  boun¬ 
daries,  may  cause  large  variations  in  fatigue  life.  Quan¬ 
titative  characterization  of  such  defects  in  a  fracture 
mechanics  basis  is  possible  by  assigning  an  initial  crack 
size  corresponding  to  the  observed  fatigue  life. 


2.  FATIGUE  TESTS  OF  FILLET  WELDED  JOINTS 


Two  cruciform  shaped  specimens  were  selected  for 
investigation  in  this  study.  Figure  1  shows  a  schematic  of 
the  two  specimens.  The  upper  specimen  is  considered  a  non¬ 
load  carrying  specimen  since  the  main  plate  is  continuous. 
Consequently,  the  welds  and  cross  plates  only  carry  a  small 
portion  of  the  load.  The  lower  specimen  has  a  discontinuous 
main  plate  and  the  full  load  must  be  carried  by  the  welds 
into  the  cross  plate. 

These  specimens  were  selected  due  to  their  ease 
of  fabrication  and  testing  at  high  stress  levels.  They  are 
representative  of  typical  joints  found  in  welded  construction 
and  have  been  used  by  numerous  investigators.  Their  use 
permitted  the  results  of  this  study  to  be  directly  compared 
with  earlier  work.  The  simple  joint  geometry  facilitated  the 
stress  analysis  and  reduced  the  number  of  geometric 
parameters . 

2 . 1  Specimen  Fabrication 

The  specimens  were  fabricated  from  5/8"  thick  ASTM 
grade  A514F  plate  material.  Table  1  summarizes  the  physical 
and  mechanical  properties.  The  specimens  were  fabricated 
from  four  sets  of  plates.  The  plates  were  cut  so  that  the 
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direction  of  rolling  of  the  main  plate  was  the  same  as  the 
direction  of  the  appl  load.  After  cutting,  the  plates 
were  sand  blasted  to  remove  loose  mill  scale.  The  surfaces 
to  be  welded  were  then  ground  lightly  to  remove  the 
remaining  tight  mill  scale.  The  plate  edges  were  left 
square. 


One-half  inch  fillet  welds  were  used  to  connect 
the  cross  plates  to  the  main  plate.  The  weld  size  was 
selected  to  provide  a  static  weld  strength  greater  than  or 
equal  to  the  capacity  of  the  main  plate.  The  welds  were 
made  in  three  passes  in  the  sequence  shown  in  Fig.  2.  They 
were  placed  in  the  horizontal  position  with  the  legs 
inclined  so  that  the  plane  of  the  electrode  travel  was  45° 
from  the  surface  of  the  two  plates. 

The  welds  were  made  by  the  automatic  gas  metal  arc 
process.  One-sixteenth  inch  AX-110  welding  wire  was  used 
with  Argon  plus  2%  oxygen  shielding  gas.  The  average  heat 
input  was  40,000  joules/in  with  an  interpass  temperature  of 
less  than  150°F.  These  welding  procedures  were  selected  to 
provide  a  weld  which  slightly  overmatched  the  strength  of 
the  base  plate. 

After  welding,  the  joints  were  sliced  into  pieces 
1  3/8  inches  wide.  The  edges  were  then  Blanchard  ground  to 


-9 


a  width  of  1.200  inches  +  0.001.  The  nominal  specimen 
dimensions  are  given  in  Fig.  1. 

Figure  3  shows  a  polished  and  etched  joint  section. 
The  three  weld  passes  are  readily  apparent  together  with  the 
heat  affected  zone  and  the  lack  of  penetration  at  the  weld 
roots.  The  variation  of  the  weld  profiles  shown  is  typical 
of  all  the  joints  tested. 

2 . 2  Testing  Procedure 

The  specimens  were  tested  in  a  220  kip  Amsler 
fatigue  testing  frame.  The  hydraulic  power  was  supplied  by 
a  220  kip  Amsler  Slow  Cycling  Unit.  The  specimens  were 
tested  under  a  constant  amplitude  cyclic  stress.  Due  to  the 
operating  characteristics  of  the  slow  cycling  unit  the 
waveform  was  a  rounded  symmetrical  saw  tooth.  The  testing 
frequency  varied  between  20  and  60  CPM.  The  lower  frequency 
was  used  for  the  highest  stress  range  specimens  and  the 
higher  frequency  lor  the  lower  stress  ranges. 

The  nominal  stress  in  the  main  plate  of  the  specimen 
was  selected  as  the  controlled  experimental  variable.  The 
minimum  stress  of  all  the  specimens  was  20  ksi.  Table  2 
gives  the  values  of  stress  range  tested  and  the  corresponding 
maximum  stress.  A  minimum  of  three  specimens  of  both  types 
of  specimens  were  tested  at  each  stress  level  to  provide  a 
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measure  of  the  variability.  The  order  of  testing  and  the 
assignment  of  specimens  to  the  stress  levels  was  randomized 
to  minimize  the  effect  of  the  uncontrolled  variables  such 
as  weld  geometry,  humidity,  and  ambient  temperature. 

Foil  strain  gages  with  a  gage  length  of  1/4"  were 
applied  to  both  sides  of  one  end  of  the  main  plate.  These 
were  located  3  inches  from  the  cross  plates.  The  gages  were 
used  to  measure  the  load  in  the  main  plate  since  the  Slow 
Cycling  Unit's  load  indication  was  not  accurate  during  cyclic 
loading.  The  strain  gage  output  was  monitored  on  an  x-y 
recorder.  The  system  allowed  the  leads  to  be  set  such  that 
stress  range  in  the  main  plate  was  within  2%  of  the  desired 
range . 


The  tests  were  run  continuously  until  either  the 
specimen  fractured  or  its  stiffness  was  reduced  suddenly  by 
a  failure  of  one  weld  causing  the  machine  to  shut  off. 

2 . 3  Failure  Modes 

The  results  of  the  fatigue  tests  are  summarized  in 
Tables  3  and  4.  The  specimen  designation  is  shown  along  with 
the  value  of  the  nominal  stress  range  and  maximum  stress  in 
the  main  plate.  The  failure  mode  of  each  specimen  and  its 
cyclic  life  is  also  shown.  Three  failure  modes  were  observed. 
All  of  the  non-load  carrying  welds  failed  from  cracks 
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extending  from  the  weld  toe  and  are  designated  as  the  toe 
mode.  Figure  4  shows  a  fracture  corresponding  to  this  mode 
of  failure.  The  crack  propagated  through  the  thickness  of 
the  main  plate  until  failure. 

Most  of  the  load  carrying  weld  specimens  failed 
from  cracks  which  propagated  from  the  initial  flaw  at  the 
weld  root.  This  discontinuity  resulted  from  the  incomplete 
penetration  of  the  weld.  These  failures  were  designated  as 
the  root  mode.  The  crack  propagated  perpendicular  to  the 
main  plate  into  the  weld  metal.  The  final  failure  occurred 
by  a  shearing  off  across  the  weld.  Figure  5  shows  a  typical 
root  failure. 

Twc  other  failure  modes  were  observed  in  the  load 
carrying  weld  joints.  Five  joints  failed  from  cracks 
extending  from  the  wold  toe  and  were  similar  to  the  failure 
shown  in  Fig.  4.  ,\  fourth  failure  mode  observed  in  throe 

load  carrying  joints  was  a  shear  failure  as  shown  in  Fig.  6. 
The  failure  plane  was  parallel  to  the  main  plate.  The 
failure  initiated  at  the  wold  toe  from  small  cracks  which 
wore  propagating  into  the  plate  until  a  few  cycles  before 
final  failure.  The  crack  growth  direction  abruptly  changed 
to  a  direction  parallel  to  the  main  plate.  The  failure 
occurring  when  the  crack  reached  the  wold's  root. 
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2.4  Test  Results 

All  test  results  are  plotted  in  Fig.  7  where  the 
life  of  the  specimen  is  plotted  as  a  function  of  the  nominal 
stress  range  in  the  main  plate.  The  non-load  carrying  joints 
are  seen  to  yield  a  longer  fatigue  life.  The  different 
failure  modes  are  indicated  by  the  symbols  in  the  legend. 

Some  of  the  welds  exhibited  marked  amounts  of  porosity  and 
are  indicated  by  vertical  lines  through  their  symbols.  They 
are  also  marked  by  astericks  in  Table  3.  The  shorter  lives 
of  load  carrying  specimens  with  porous  welds  is  seen  to  be 
substantial  at  the  higher  stress  ranges. 

2.4.1  Effect  of  Weld  Angle  on  the  Cracking  Behavior  of  the 
Non-Load  Carrying  Weld  Joints. 

All  non-load  carrying  weld  joints  failed  from 
cracks  starting  at  the  weld  toe.  Each  joint  had  four  weld 
toes  that  provided  potential  failure  sites.  Photoelastic 
studies  have  shown  that  the  stress  at  the  weld  toe  is  a 
function  of  the  weld  angle.  '  The  influence  of  the  weld 
angle,  0,  on  the  failure  location  and  specimen  life  was 
examined  to  determine  its  affect.  The  weld  angle  was  defined 
by  the  weld  profile  as  shown  in  Fig.  8.  The  measurements  of 
the  weld  height  (II)  and  weld  width  (W)  were  made  with  a 
machinist's  scale  to  the  nearest  hundredth  of  an  inch. 
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Figure  9  provides  a  histogram  of  the  measured  weld 
angles.  The  mean  angle  of  all  four  welds  of  each  joint  is 
seen  to  be  less  than  the  mean  weld  anqle  where  failure 
occurred.  This  difference'  suggests  that  the  location  of  fail 
ure  for  a  particular  specimen  is  related  to  the  w'eld  ancle. 

The  distribution  of  the  ratio  of  each  weld  angle 
in  a  specimen  to  the  failure  weld  angle  is  given  in  Fig.  10. 
Only  6  of  the  weld  angles  were  greater  than  the  weld  angle 
at  the  weld  where  failure  occurred.  It  is  apparent  that  the 
larger  the  angle  of  a  weld  the  more  likely  the  failure  will 
occur  at  that  weld.  Ninety-four  percent  of  the  joints  in 
this  test  series  failed  at  welds  with  the  largest  angle. 

The  fact  that  the  failure  location  for  a  joint  is 
dependent  upon  the  angle'  of  the  welds  sugacsts  that  the 
fatigue  life  is  influenced  by  the  weld  angle.  Figure  11  was 
constructed  to  examine  this  hypothesis.  The  ordinate  defines 
the  ratio  of  the  weld  angle  at  the  failure  location  to  the 
mean  of  the  weld  angles  where  failure  occurred.  The  abscissa 
was  taken  as  the  ratio  of  the  life  of  the  specimen  to  the' 
average  life  for  that  stress  range.  Although  substantial 
scatter  is  apparent,  a  trend  between  weld  angle  and  life  is 
apparent.  The  solid  lino  is  a  least  square  fit  to  the 
plotted  data.  The  dashed  line  represents  the  mean  line  less 
the  standard  error  of  estimate.  The  life  of  a  particular 
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specimen  is  therefore  a  function  of  the  weld  angle  as  well 
as  the  applied  stress. 

The  scatter  in  Fig.  11  seems  reasonable  since  the 
weld  angle  is  only  a  gross  indication  of  the  weld  profile. 
Small  irregularities  at  the  toe  along  the  length  of  the  weld 
such  as  undercut  also  affect  the  stress  condition  at  the 
weld  toe  and  the  initial  crack  condition  and  hence  the 
fatigue  life.  The  weld  angle  defined  in  this  study  provides 
a  rough  estimate  of  the  variation  in  the  stress  at  the  weld 
t<_>e  when  the  weld  profile  changes. 

The  influence  of  the  weld  angle  upon  the  fracture 
path  of  the  toe  failures  was  also  investigated.  The  angle, 
a,  which  the  crack  made  with  the  vertical  was  measured  as 
shown  in  Fig.  8.  The  measurements  were  made  with  a  pro¬ 
tractor  to  an  estimated  accuracy  of  +1°.  An  attempt  was 
made  to  correlate  the  crack  angle  with  the  weld  angle  since 
the  state  of  stress  at  the  toe  is  dependent  upon  the  weld 
angle.  The  crack  plane  was  expected  to  be  perpendicular  to 
the  direction  of  the  principal  stress.  Figure  12  shows  a 
plot  of  the  weld  angle  and  crack  angle.  No  correlation  is 
apparent.  An  average  angle  of  13°  seems  to  be  a  reasonable 
estimate  for  all  weld  angles. 
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2.4.2  Fracture  Path  of  Toe  Failure 

Failures  from  cracks  at  the  weld  toe  occurred  in 
both  the  load  carrying  and  non-load  carrying  joints.  The 
origin  of  the  crack  was  the  same  in  both  joints.  The  cracks 
were  observed  to  start  from  the  intersection  of  the  fusion 
line  with  the  weld  surface.  Figures  13  and  14  are  photo 
micrographs  of  typical  toe  cracks.  Crack  growth  at  the 
fusion  line  is  visible  in  both  figures.  Figure  13  is  a 
fractured  specimen  that  was  tested  at  a  stress  range  30  ksi. 
The  crack  grew  away  from  the  weld  and  the  fusion  line  into 
the  base  metal  of  the  main  plate.  The  crack  shown  in  Fig. 

14  is  from  a  specimen  tested  at  a  stress  range  of  80  ksi. 

The  crack  followed  the  fusion  line  and  wash  up  zone  into 
the  weld  metal.  The  fracture  path  followed  inclusions  that 
were  observed  in  this  area.  The  two  fracture  paths  appear 
to  be  stress  dependent.  Specimens  tested  at  higher  stress 
ranges  showed  a  distinct  tendency  for  the  crack  to  follow 
the  fusion  line.  In  the  specimens  tested  at  the  lower  stress 
levels  the  crack  grew  immediately  into  the  base  plate.  At 
the  higher  stress  levels  the  crack  followed  the  fusion  line 
for  a  considerable  distance  before  changing  direction  and 
growing  in  the  plate.  Figure  15  shows  the  specimen  which 
exhibited  the  longest  fusion  line  growth  of  all  the  specimens 
tested.  In  most  of  the  specimens  where  crack  growth  was 
observed  along  the  fusion  line,  the  crack  length  along  the 
fusion  line  was  of  the  order  of  0.05  inches.  Figure  16  shows 
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an  enlargement  of  the  area  where  the  crack  deviated  from  the 
fusion  line  into  the  plate. 

The  tendency  for  the  crack  to  follow  the  fusion 
line  at  high  stress  levels  seems  to  be  related  to  the  inclu¬ 
sion  content  of  the  heat  affected  zone.  Figure  17  shows  an 
unetched  section  of  the  weld  metal  and  heat  affected  zone. 

A  large  number  of  inclusions  are  visible  in  the  heat  affect¬ 
ed  zone.  Several  elongated  inclusions  are  parallel  to  the 
fusion  line.  An  examination  of  the  heat  affected  zone  at  a 
higher  magnification  can  be  made  from  the  etched  specimen 
shown  in  Fig.  18.  The  inclusions  at  the  grain  boundaries 
are  easily  seen.  Heiser  has  shown  that  in  banded  steels  the 

inclusions  play  an  important  role  in  fatigue  crack  propaga- 
(13) 

tion.  At  high  stress  levels  he  found  that  the  inclusion 

content  had  a  significant  effect  upon  fatigue  crack  growth 
rates  and  fracture  paths.  The  crack  tended  to  follow  inclu¬ 
sions  and  grew  intergranularly  between  the  included  grain 
boundary  particles.  He  suggested  that  this  behavior  was  due 
to  the  inability  of  the  included  particles  to  adjust  to  the 
stress  conditions  causing  cracking  around  the  particle. 

This  cracking  left  a  void  which  the  crack  propagated  along. 
At  lower  stress  conditions  cracking  at  included  particles 


was  not  observed. 
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The  high  inclusion  content  of  the  area  adjacent 
to  the  fusion  line  seems  to  be  responsible  for  the  crack 
following  that  path.  Figure  19  shows  a  crack  which  formed 
at  the  toe  and  propagated  in  the  HAZ  following  the  fusion 
line.  The  crack  branching  at  inclusions  is  visually  evident. 
Figure  20  is  a  higher  magnification  of  the  lower  part  of 
Fig.  19. 


This  study  has  shown  that  the  length  of  the  frac¬ 
ture  path  along  the  fusion  line  is  dependent  upon  the  stress 
level.  The  higher  the  stress  level  the  longer  the  crack 
followed  this  path  before  it  grew  through  the  thickness  of 
the  main  plate.  This  behavior  was  attributed  to  inclusion 
cracking  at  high  stress  levels.  The  larger  inclusion  con¬ 
tent  of  the  HAZ  at  the  fusion  line  caused  the  crack  to 
follow  this  path  at  higher  stress  levels.  All  toe  cracks 
initiated  at  the  intersection  of  the  fusion  line  with  the 
plate  surface. 

2.4.3  Fracture  Path  of  Root  Failures 

Root  failures  only  occurred  in  the  load  carrying 
weld  joint.  No  visible  evidence  of  crack  growth  was  found 
at  the  weld  roots  of  the  non-load  carrying  welds. 

Figures  21  and  22  show  the  region  at  the  weld  root  for  two 
untested  load  carrying  joints.  The  largo  separation 
at  the  riaht  of  the  photos  is  due  to  the  fit  up  of  the 
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plates.  The  sharpness  of  this  separation  in  Fig.  21  is 
seen  to  increase  as  it  approaches  the  weld  metal.  Both 
figures  indicate  an  extremely  sharp  crack  immediately 
adjacent  to  weld  metal.  Cracks  which  caused  the  root 
failures  propagated  from  these  initial  cracks. 

Figures  23  and  24  show  typical  fracture  paths  of 
the  root  crack  in  the  weld  metal.  The  crack  path  has  fol¬ 
lowed  the  columar  grain  boundaries  in  Fig.  23.  This  ten¬ 
dency  is  also  evident  in  Fig.  24  but  to  a  lesser  degree. 

The  path  of  the  crack  is  nearly  perpendicular  to  the  main 
plate  in  both  figures.  A  tendency  for  the  crack  to  follow 
the  grain  boundaries  in  the  weld  metal  was  observed  at  all 
the  root  cracks  examined.  However  the  degree  to  which  this 
influenced  the  macroscopic  fracture  path  was  found  to  be 
negligible.  The  crack  would  follow  the  grain  boundaries 
when  they  were  favorably  oriented  with  respect  to  the 
macroscopic  fracture  path.  If  the  grains  of  the  weld  were 
oriented  perpendicular  to  the  main  plate  at  the  weld  root, 
as  shown  in  Fig.  23,  the  fracture  path  was  along  the  grain 
boundaries.  When  the  grains  were  not  perpendicular  to  the 
oriented  plate,  the  crack  would  grow  transgr.mularly  across 
the  grains  as  shown  in  Fig.  24. 

These  cracks  grew  vertically  into  the  weld  metal 
until  failure  occurred  by  a  sudden  tear  as  shown  in  Fig.  6. 


Estimates  of  the  stress  on  the  final  failure  plane  were 
made  and  found  to  be  smaller  than  the  ultimate  strength  of 
the  weld  metal.  This  fact  suggests  that  weld  metal  may 
have  experienced  fatigue  damage  during  testing.  This  caused 
the  fracture  to  occur  at  a  lower  stress  than  its  ultimate 
capacity.  No  evidence  of  such  damage  was  found ,  although 
as  shown  in  Fig.  15  the  failure  across  this  plane  may  have 
been  due  to  rapid  crack  growth. 

Figure  25  shows  the  tip  of  a  root  crack  at  failure 
after  the  crack  has  changed  direction.  The  large  open  crack 
is  perpendicular  to  the  main  plate  while  the  smaller  branch 
is  almost  parallel  to  the  main  plate.  The  final  failure, 
therefore,  was  due  to  rapid  fatigue  crack  extension.  The 
change  in  the  fracture  path  was  attributed  to  a  change  in 
the  principal  strops  direction  caused  by  the  bending  moment 
introduced  in  the  weld  when  the  root  crack  is  large. 

The  porosity  of  the  weld  also  had  a  dramatic  effect 
upon  the  fatigue  lives  of  specimens  which  failed  from  root 
cracking.  Figure  7  shows  that  at  high  stress  levels 
porosity  in  the  weld  caused  failure  to  occur  at  much  shorter 
lives.  This  effect  is  believed  to  be  due  to  crack  growth 
at  the  porosity.  Crack  growth  was  evident  at  the  porosity 
in  all  specimens  with  porous  welds.  These  cracks  weakened 
the  weld  and  decreased  its  faituge  life.  Figure  26  shows 
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the  failure  of  a  joint  with  a  porous  weld.  The  failure  oc¬ 
curred  by  the  rupture  of  the  weld  along  the  ligaments  between 
the  extending  cracks.  At  lower  stress  levels  the  load  on 
these  ligaments  was  less  and  failure  was  similar  to  the  root 
failures  of  welds  without  porosity. 

2.4.4  Analysis  of  Shear  Type  Failures  of  the  Load  Carrying 
Weld  Joint 

Three  load  carrying  welds  exhibited  a  shear  type 
failure  mode  as  shown  in  Fig.  6.  Examination  of  the  welds 
revealed  that  the  leg  of  the  weld  (W  in  Fig.  1)  was  much 
smaller  than  the  average.  This  fact  along  with  the  mode  of 
failure  suggested  that  the  failure  of  these  specimens  was 
due  to  high  shear  stresses  on  the  base  of  the  weld.  Figure 
27  shows  the  lives  of  the  load  carrying  welds  plotted 
against  the  shear  stress  range  of  the  weld  with  the  smallest 
leg.  The  stress  was  calculated  by  dividing  the  total  load 
range  per  unit  width  by  the  smallest  weld  leg.  The  line 
shown  in  the  figure  is  a  least  square  fit  of  the  shear 
failure  data. 

The  toe  failures  at  high  stress  levels  are  the  only 
points  above  the  least  square  line.  These  failures  occurred 
by  crack  extension  into  the  main  plate  after  following  the 
fusion  line  for  a  portion  of  their  life.  The  shear  failures 
were  similar  to  the  toe  failure  except  that  the  crack  did 
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not  propagate  into  the  main  plate.  The  failure  occurred  by 
cracking  along  the  fusion  line.  At  lower  stress  levels  less 
cracking  was  evident  along  the  fusion  line.  The  line  in 
Fig.  27  provides  a  reasonable  estimate  of  the  fatigue  life 
of  the  weld  in  shear.  At  high  stress  levels,  the  crack 
propagated  along  the  fusion  line  due  to  the  inclusion 
cracking  until  the  stress  conditions  favored  through  the 
plate  thickness  growth  (toe  failures) . 

If  the  inclusion  orientation  and  spacing  are 
favorable  the  crack  may  continue  along  the  fusion  line 
resulting  in  a  shear  mode  failure.  The  tendency  for 
inclusion  cracking  is  less  at  lower  stress  levels  and  hence 
the  specimens  were  not  as  critical  in  shear.  Other  failure 
modes  were  more  critical  as  is  apparent  in  Fig.  27. 
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3.  FRACTURE  MECHANICS  ANALYSIS  OF  FATIGUE  BEHAVIOR 


Fracture  mechanics,  developed  from  the  study  of 
brittle  fracture,  has  found  increasing  use  in  the  study  of 
structural  fatigue. '  '  1  Fracture  mechanics  is  used  to 

relate  the  materials  fracture  behavior  by  considering  the 
energy  balance  of  the  energy  released  in  fracture  to  the 
materials  ability  to  absorb  this  energy.  The  use  of 
fracture  mechanics  in  structural  fatigue  is  based  on  the 
assumption  that  the  fatigue  crack  growth  rate  is  a  function 
of  the  crack  tip  stress  field. 

3. 1  Crack  Tip  Stress  Field  Equations 

The  stress  field  in  the  region  immediately 
adjacent  to  the  crack  tip  for  a  two  dimensional  body  is  of 
the  form: 

CT  cos9/2[l-sin9/2sin38/2]  -sin9/2[2+cos9/2cos39/2] 

X 

1  KI 

Cy  *=  c03®/2[l+sin9/2sin38/2}  sin9/2cos9/2cos39/2  (1) 

!  ,  ,  i  Kn 

T  |  sin8/2cosO/2cos3B/2  sin6/2cos8/2cos38/2 

xy  l —  j 

where  the  coordinates  "r"  and  6  are  shown  in  Fig.  28.  The 

parameters  and  Kj  ,  called  the  crack  tip  stress-intensity 

factors,  are  two  independent  parameters  which  are  functions 

of  the  body's  geometry,  crack  length  and  linear  functions  of 


r 
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the  applied  stress.  represents  the  contribution  of  a 
symmetric  and  an  antisymmetric  stress  distribution  with 
respect  to  the  crack  plane  on  the  crack  tip  stress  field. 
The  stress  intensity  factors  are  independent  of  the 
coordinates  "r"  and  0  hence,  represent  the  intensity  of  the 
stress  field. 


The  stress-intensity  factors  are  the  leading  terms 
of  an  expansion  of  the  infinite  series  representing  the 
total  stress  distribution  of  the  body.  The  stress  field 
adjacent  to  the  crack  is  dominated  by  a  square  root  singu¬ 
larity  in  "r"  at  the  crack  tip.  The  displacement  field  at 
the  crack  tip  is  also  dependent  on  the  stress  intensity 
factor.  Since  the  fatigue  process  takes  place  at  the  crack 
tip,  it  is  reasonable  to  assume  that  it  must  be  a  function 

of  the  crack  tip  stress  field  and  consequently  the  stress- 

(35' 

intensity  factors.  ; 


3. 2  Relationship  Between  the  Stress-Intensity  Factor  and 
Fatigue  Crack  Growth  Rate 


Paris  proposed  the  following  relationship  between 
the  stress  intensity  range  (AK)  and  the  crack  growth  rate 
(da/dN)  (36) 


da/dN  =  CAKn 


(2) 


This  relationship  is  a  semi-empirical  attempt  to  fit  experi¬ 
mental  data  with  the  stress  intensity  factor.  Numerous 
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crack  growth  studies  on  a  variety  of  materials  have 
shown  that  Eq.  2  provides  a  reasonable  estimate  of  a 
material's  crack  growth  rate  response. (33,37,38)  other 
equations  similar  to  Eq.  2  have  been  proposed.  Some  of 
these  equations  have  been  based  on  work  hardening  and 
fatigue  damage  considerations  at  the  crack  tip  region.  Other 
equations  have  been  based  on  consideration  of  the  crack  tip 
displacements.  Reference  39  gives  an  up-to-date  review  of 
these  different  models.  Most  models  reduce  to  the  general 
form  of  Eq.  2  for  a  particular  material  and  environment. 

3. 3  Analysis  of  Fatigue  Life  Behavior  Using  Fracture 
Mechanics 

The  power  law,  (Eq.  2) ,  relating  the  fatigue  crack 
growth  rate  to  a  power  of  the  range  of  stress-intensity 
provides  the  necessary  relationship  to  estimate  the  fatigue 
life  of  a  component.  The  range  of  the  stress-intensity 
factor  AK  may,  in  general,  can  be  expressed  as 

AK  =  SR  g (a,w' )  (3) 

where  S_  =  nominal  stress  range 

K 

a  =  crack  length 

w'  =  a  set  of  geometric  constants  for 
the  body. 

When  Eq.  3  is  substituted  into  Eq.  2  it  yields 

g  (a ,w' ) ”nda  =  C  SRndN  (4) 

Integratinn  Eq.  4  between  the  limits  of  the  initial  crack 
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length  and  the  final  crack  length  a^,  defines  the 
component's  fatigue  life  N  as 

."3  f 

N  =  )a_  g (a,w' ) "nda/C  SRn  (5) 

Equation  5  provides  a  means  to  estimate  the 
fatigue  life  of  a  component  if  the  value  of  SR,  a^,  a^,  C, 
n,  and  the  function  g  are  known.  It  is  also  assumed  that  the 
fatigue  behavior  of  the  component  is  totally  a  function  of 
crack  growth.  Crack  initiation  is  not  considered. 

It  is  also  interesting  to  note  that  the  values  of 
C  and  n  can  be  estimated  from  Eq.  5.  If  the  values  of  a^ 
af,  N,  and  SR  are  known  from  a  series  of  fatigue  tests,  and 
the  function  "g"  is  known  for  the  geometry  of  the  specimen, 
then  the  values  of  C  and  n  can  be  determined  by  fitting  the 
results  of  the  fatigue  tests  to  Eq.  5. 
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This  provides  a  linear  relationship  of  the  form 

log  N  =  log  B  -  m  log  SR  (8) 

or  N  =  B  SR'm  (9) 

The  constant  B  can  be  expressed  as 

B  =  X/C  (10) 

when  the  exponents  "n"  and  "m"  are  equated.  Equation  6 
provides  a  means  of  including  the  effect  of  component 
geometry  and  initial  flaw  size,  a^,  by  the  term  n/X  .  The 
consequences  of  this  equation  are  explored  and  analyzed  for 
fillet  welded  joints  in  this  dissertation. 
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4.  ESTIMATION  OF  STRAIN  ENERGY  RELEASE  RATES 
USING  THE  FINITE  ELEMENT  TECHNIQUE 


Irwin  has  shown  that  the  stress-intensity  factor, 

K,  is  related  to  the  strain  energy  release  rate,  G.  For 

( 34) 

a  plane  stress  analysis  this  can  be  expressed  as 

G  =  K:2  /  E  +  KXI 2  /  E  (11) 

The  value  of  the  strain  energy  release  is  defined  as  the 

energy  rate  provided  by  a  body  containing  an  extending 

crack  available  for  the  crack-extension  process.  In  other 

words  the  energy  per  unit  of  new  crack  area  generated.  This 

energy  rate  can  be  related  to  the  change  in  compliance  of 

( 34 ) 

the  body  for  a  load  P  as 

G  =  P2/2  dC/dA  (12) 

where  dC/dA  is  the  rate  of  change  of  compliance  with  crack 
area  where  the  compliance  is  defined  as 

C  =  A/P  (13) 

4 . 1  Finite  Element  Analysis 

The  finite  clement  technique  has  found  wide  spread 
use  in  the  solution  of  boundary  value  problems.  The 
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method  is  of  a  general  nature  and  can  be  applied  to  any 
boundary  value  problem. 

The  method  of  analysis  can  be  divided  into  two 

stages : 

1.  The  discretizing  of  the  body  to  be  solved  into 
finite  elements,  and 

2.  Minimization  of  a  selected  functional  in  the 
body . 

The  first  stage  in  the  analysis  is  done  to  make  the  second 
stage  more  tractable.  By  dividing  the  body  into  discrete 
elements  the  functional  to  be  minimized  is  represented 
continuously  in  each  element.  The  result  is  a  piecewise 
continuous  representation  of  the  functional  over  the  body. 
The  discrete  elements  are  connected  at  nodes  with  the  fields 
in  each  element  written  as  functions  of  nodal  parameters. 

The  boundary  value  problem  is  therefore  reduced  to  the 
solution  of  a  set  of  simultaneous  equations  in  terms  of  the 
nodal  parameters  connecting  the  finite  elements.  These 
sets  of  equations  are  of  the  form  of 

if)  =  [A  J  16}  (14) 

where  "f"  is  a  vector  of  nodal  forces  and  "6"  the  nodal 
displacements.  The  matrix  "A"  is  the  stiffness  matrix  of 
the  discretized  body  relatinq  these  parameters.  The 
accuracy  of  a  finite  element  solution  depends  on  how  well 
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the  piecewise  representation  of  the  functional  from  the 
analysis  approximate  its  true  form. 

4 . 2  Finite  Element  Compliance  Analysis  of  Cruciform 
Weldments  ~ 

The  compliance  of  the  cruciform  joint  was  deter¬ 
mined  from  a  finite  element  analysis.  The  body  was 
analyzed  as  a  two  dimensional  plane  stress  problem.  Figure 
29  shows  a  schematic  of  the  idealized  shape  of  the  two 
joints  analyzed.  The  welds  profiles  were  idealized  as 
triangular  shapes.  The  weldments  were  also  assumed  to  be 
completely  symmetric  with  respect  to  the  planes  shown  as 
dashed  lines  in  Fig.  29.  The  length  of  the  main  plate  from 
the  weld  toe  was  taken  as  three  times  the  thickness  of  the 
main  plate.  This  located  the  loaded  nodes  far  enough  away 
from  the  weld  toe  and/or  toe  cracks  so  that  the  state  of 
stress  at  those  points  would  be  uniform  axial  tension. 
Symmetry  of  the  joint  permitted  a  quarter  section  of  the 
total  joint  to  be  used. 

4.2.1.  Element  and  Mesh  Selection 

The  selection  of  the  type  of  elements  to  be  used 
and  the  mesh  or  discretizing  of  the  body  determines  the 
accuracy  of  the  finite  element  solution.  The  more  sophis¬ 
ticated  the  element  the  larger  the  mesh  can  be  for  the  same 
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accuracy.  A  more  sophisticated  element  contains  a  higher 
order  field  and  consequently  can  represent  the  functional 
over  a  larger  area. 

The  crack  tip  stress  field  equations  given  in 
Chapter  3  show  that  at  the  crack  tip  a  stress  singularity  of 

the  form  l//r  exists.  Consequently  a  displacement  field 

x 

of  the  form  r2  exists  in  the  region  of  the  crack  tip.  Due 
to  the  geometry  of  the  weld  joint,  a  complex  stress  distri¬ 
bution  exists  in  the  region  of  the  weld  and  at  the  weld's 
toe  and  root . 

The  experimental  study  and  the  work  of  others  have 
shown  that  cracking  can  occur  at  two  locations  in  load 
carrying  weld  joints.  Cracks  may  occur  at  the  root  of  the 
weld  or  at  the  weld  toe  depending  upon  the  geometry  of  the 
joint  and  the  penetration  of  the  weld.  In  non-load 
carrying  weld  specimens,  cracking  nearly  always  occurs  at 
the  weld  toe.  Hence,  a  finite  element  analysis  of  load 
carrying  joints  must  provide  for  too  and  root  cracks.  Only 
toe  cracks  need  to  he  considered  for  the  non-load  carrying 
weld  joints. 

Due  to  the  complexity  of  t  lit'  stress  distribution 
and  the  large  number  of  geometries  and  crack  lengths  to  be 
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analyzed,  an  element  was  selected  that  would  minimize  the 
time  and  cost  of  the  analysis.  The  Constant  Strain 
Triangle  (CST)  element  is  the  simplest  finite  element 
available  and  was  selected.  Its  chief  advantages  are  that 
its  element  properties  are  easily  and  rapidly  calculated 
and  its  triangular  shape  facilitated  the  matching  of 
geometric  boundaries. 

The  displacement  field  within  the  element  is 
linear  so  consequently  the  strain  and  stress  field  are 
constant.  Each  element  provides  a  constant  strain 
energy  density.  Since  the  strain  energy  was  selected  as 
the  functional  to  be  minimized  (as  is  normally  done  in  the 
solution  of  boundary  value  problems  whose  solution  is 
formulated  from  displacement  fields)  it  is  approximated 
stepwise  by  each  element  over  the  body. 

The  element  mesh  and  size  in  the  regions  of  high 
displacement  gradients  must  be  refined  to  insure  that  the 
true  displacement  field  or  strain  energy  density  is 
adequately  represented.  figure  10  shows  the  mesh  used  to 
evaluate  thM  load  carrying  weld  with  a  root  crack .  The 
element  size  near  the*  crack  tip  and  in  the  region  of  the 
weld  toe  is  seen  to  be  quite  small.  The  elements  within 
the  weld  were  also  made  small  to  accommodate  the  complex 
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state  of  stress  which  existed.  Figures  31  and  32  show  the 
meshes  used  for  the  analysis  of  toe  cracks.  Again  the 
refinement  in  the  mesh  at  the  crack  tip  is  evident. 

The  crack  in  each  of  the  joints  was  modeled  by 
two  lines  of  nodal  points  immediately  adjacent  to  each  other 
but  not  connected.  Displacement  boundary  conditions  were 
applied  along  the  planes  of  symmetry.  The  nodes  at  the 
bottom  of  the  figures  at  the  mid-thickness  of  the  main 
plate  were  fixed  in  the  vertical  direction.  The  nodes  at 
the  mid-thickness  of  the  cross  plate  were  restrained  in  the 
horizontal  direction.  The  nodes  at  the  end  of  the  main 
plate  were  loaded  by  concentrated  tension  loads  equivalent 
to  a  uniform  stress  distribution. 

The  small  elements  used  in  the  mesh  near  the  crack 

tip,  although  much  smaller  than  normally  used  in  finite 

element  analysis  of  bodies,  arc  still  larger  than  those 

used  by  previous  investigators  of  cracked  bodies .  ^  °  '  **  *  ' 4  ^  ^ 

This  was  not  believed  to  invalidate  the  results  since  the 

values  of  the  stress  intensity  factors  were  determined  from 

a  compliance  analysis.  Mowbray  has  recently  demonstrated 

that  this  method  gives  quite  accurate  results  with  rather 
(43) 


coarse  meshes. 
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The  reason  for  this  accuracy  is  due  to  the  fact 

that  although  the  finite  element  solution  underestimates  the 

strain  energy  near  the  crack  tip,  the  compliance  method  of 

analysis  examines  the  strain  energy  in  the  entire  body. 

Since  the  crack  tip  singularity  only  dominates  the  region 

close  to  the  tip  the  overall  estimate  of  the  strain  energy 

is  good.  The  strain  energy  at  the  crack  tip  region  is  only 

a  small  percentage  of  the  total  strain  energy  in  the  body. 

Consequently  errors  in  its  estimation  do  not  have  much  effect 

on  the  total  strain  energy  calculation.  Methods  suggested 

by  other  investigators  focus  on  the  local  behavior  at  the 

crack  tip  and  consequently  require  extremely  small  elements 

to  accurately  estimate  the  stress  and  displacement  behavior 
(41  42  43) 

in  this  region.  '  '  High  order  elements  which  include 

the  stress  singularity  at  the  tip  do  not  require  the 

refinement  in  mesh  size,  however  displacement  discontinuities 

arise  between  these  elements  and  adjacent  ones  which  arc 

(44) 

also  a  source  of  error. 

The  effect  of  the  error  in  the  estimation  of  the 
strain  energy  near  the  crack  tip  on  the  analysis  can  also  be 
evaluated  by  comparison  of  compliance  tests  of  actual 
specimens  with  analytical  solutions.  Srawley  found  the  error 

(45) 

to  be  small  for  a  single  edge-notch  specimen.  The 

compliance  was  measured  on  a  specimen  where  the  crack  tip 
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radius  was  1/32  inch.  The  stress  at  the  crack  tip  is 
proportional  to  the  ratio  of  the  crack  length  to  the  crack 
tip  radius  with  the  value  of  the  stress  intensity  factor 
being  defined  as  the  crack  tip  radius  approaches  zero. 
Consequently  the  stress  field  in  the  vicinity  of  the  crack 
tip  was  less  than  that  defined  by  equation  1.  This  is 
similar  to  the  effect  of  the  underest iuiate  of  this  stress 
field  by  the  finite  element  solution.  Therefore, 
the  assumption  that  the  error  in  estimating  this  energy  will 
not  influence  the  accuracy  of  the  compliance  analysis  seems 
justified. 


The  second  reason  given  by  Mowbray  for  the  accuracy 
of  the  method  is  the  fact  that  the  values  of  G  and  K  are 
functions  of  the  change  in  compliance  with  crack  length. 
Consequently,  if  the  error  in  compliance  is  independent  of 
crack  length,  the  slope  of  the  compliance  versus  crack  length 
curve  should  be  close  to  the  true  curve.  The  independence 
of  the  error  with  crack  length  should  hold  if  the  finite 
element  mesh  is  similar  for  each  crack  length  analyzed  as  it 
was  in  this  study. 

4.2.2  Method  of  Solution  of  Finite  Element  Equations 

The  method  of  solution  of  the  equations  resulting 
from  the  finite  element  formulation  was  selected  to  provide 
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rapid  solution  time.  The  quantities  required  for  a 
solution  by  a  compliance  analysis  are  the  displacements  of 
the  loaded  nodes.  The  loads  at  these  nodes  were  based  on 
the  assumption  that  the  state  of  stress  was  uniform.  The 
applied  stress  does  not  influence  the  resulting  load- 
displacement  relationship  so  for  convenience  in  checking 
the  results  it  was  taken  as  100  ksi  for  all  geometries. 

The  load-displacement  relationship  was  found  by  treating  the 
joint  as  a  series  of  substructures.  The  load  displacement 
relationship  was  determined  for  each  substructure  and  then 
combined  to  give  the  load-displacement  relationship  for  the 
joint.  Appendix  A  gives  the  details  of  the  analysis.  This 
method  reduces  the  size  of  the  matrix  inversions  required 
and  eliminates  most  of  the  off  diagonal  zero  terms.  For 
each  crack  length  analyzed,  only  the  substructures  near  the 
crack  are  subjected  to  mesh  changes,  consequently  only 
those  substructure  stiffness  matrices  need  be  recomputed  for 
each  crack  length. 

Each  analysis  provided  the  displacements  of  the 
loaded  nodes.  A  set  of  displacements  was  determined  for 
each  joint  noometry  corresponding  to  a  set  of  crack  lengths. 
The  compliance  of  a  particular  geometry  and  crack  length  was 
then  determined  by  dividing  the  load,  P,  in  the  main  plate 
by  the  x-disp lacement  of  the  loaded  node  at  the  plane  of 
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symmetry  of  the  main  plate.  The  x-displacement  was 

doubled  since  it  represents  the  total  extension  of  the 
symmetric  section. 


4 .  3  Evaluation  of  the  Stress  Intensity  Factor 


The  set  of  compliances  and  crack  lengths  were 
fitted  by  least  squares  to  a  polynomial  of  the  following 
form 


EBC  =  A  +  £  A  ,  (a/w')n 

n=2 


(15) 


where 

C  =  P/A  =  specimen's  compliance, 

B  =  specimen's  width, 

E  =  Young's  Modulus, 
a  =  crack  length 

and  w'=  a  geometric  constant  with  dimensions  of 
length.  The  elastic  Modulus  and  the  specimen's  thickness 
were  constant  and  taken  equal  to  30,000  ksi  and  unity 
respectively.  The  first  order  term  of  the  polynomial  was 
not  included  since  the  stress  intensity  factor 
must  be  zero  when  the  crack  length  is  zero. 


Equation  12  can  be  rewritten  as 


G 


8F.B2w ' 


d  (EBC) 


(16) 


d (a/w ' ) 
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where  the  area  of  the  new  crack  formed  is 

dA  =  4Bw '  d ( a/w ' )  (17) 

The  constant  four  in  Eq.  17  is  a  consequence  of  the  symmetry 
used  in  the  finite  element  analysis.  Each  specimen  contain¬ 
ed  four  cracks. 

Substituting  Eq.  15  into  Eq.  16  yields 

G  =  P2/8EB2w '  2  n  A  ,  (a/w')n_1  (18) 

n=2  n~1 

where  G  is  the  total  strain  energy  release  rate.  The  value 

of  the  stress  intensity  factor  was  evaluated  using  Eq.  18  as 

K  =  [EG]2  =  P/2B/2wr  2  [nA  .  (a/w')  L]  (19) 

n=2 

where  K  is  an  equivalent  stress  intensity  factor  equal  to 

2  -- 

Kj [1+ (Kjj/Kj)  ]2.  The  displacements  of  the  nodes  along  the 
crack  surface  in  the  finite  element  analysis  revealed  that 
the  ratio  of  to  Kj  was  less  than  10  percent.  Therefore 
the  value  K  calculated  by  Eq.  19  is  approximately  equal  to  Kj. 
The  stress  intensity  factors  calculated  from  the  finite  ele¬ 
ment  analysis  are  designated  as  K  in  this  dissertation. 

4.3.1  Stress  Intensity  Factors  for  Root  Cracks  of  Load 
Carrying  Fillet  Welds 

Thirty-two  different  joint  geometries  were 
analyzed  to  determine  the  effect  of  the  joint  geometry  upon 
the  stress  intensity  factor.  The  compliance  of  each  geometry 
was  determined  for  sixteen  different  crack  lengths  and  the 
uncracked  joint  (full  root  penetration).  Nine  of  the  joint 
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geometries  simulated  the  average  geometry  of  the  fatigue 
test  specimens.  The  remaining  geometric  configurations 
correspond  to  specimens  tested  in  previous  investigations. 
Table  5  summarizes  the  dimensions  of  the  joints  analyzed 
and  their  source.  The  dimensions  are  defined  in  Fig.  1. 

The  analysis  of  the  joints  simulating  the  fatigue 
test  specimens  included  a  variation  of  the  weld  angle,  0. 

The  weld  size  was  defined  by  the  leg  dimensions  W  and  H. 

These  parameters  were  determined  by  maintaining  the  weld 
cross  sectional  area  constant  and  varying  the  weld  angle. 

The  weld  area  was  equated  to  the  average  area  of  the 
specimens.  One  specimen  was  analyzed  with  an  oversize  weld 
with  equal  legs.  It  is  identified  with  an  asterisk  in  Table 
5.  The  weld  angle  was  assumed  equal  to  45  degrees  for  the 
earlier  studies  since  its  value  was  not  reported. 

The  results  for  three  test  specimens  are  plotted 

in  Fig.  33.  The  quantity  EBC  is  shown  as  a  function  of 

a/w'  for  specimens  with  weld  .ingles  of  30,  40  and  45  degrees. 

The  parameter  u'  was  defined  as 

w'  =  T  /2  +  H  (20) 

P 

for  all  the  geometries  analyzed.  A  fifth  order  polynomial 
of  the  form  of  l’q.  15  was  used  to  fit  the  compliance  (EBC) 
as  a  function  of  a'w'  for  all  geometries. 


This  order 
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polynomial  was  found  to  provide  the  minimum  standard  error 
of  estimate.  The  standard  error  of  estimate  of  the  fitted 
polynomial  was  less  than  0.004  for  all  geometries  which 
corresponds  to  an  error  of  less  than  one- tenth  of  a  percent. 

The  values  of  a/w*  at  which  the  compliance  was 
determined  were  the  same  for  all  geometries.  The  crack 
lengths  included  an  uncracked  joint  (a/w'=0)  and  increments 
of  a/w'  between  0.115  and  0.692.  These  crack  lengths  were 
believed  to  be  in  the  most  useful  range  based  on  visual 
examination  of  the  test  specimens.  The  fitted  polynomial 
provided  a  means  of  interpolating  the  value  of  the  stress 
intensity  factor  for  crack  lengths  smaller  than  w'  x  0.115. 
Accurate  calculation  of  the  compliance  for  cracks  of  this 
size  would  have  required  extremely  small  elements  using  the 
finite  element  technique.  Stress  intensity  factors  were 
determined  from  Eq.  19.  The  fitted  polynomial  permitted  the 
stress  intensity  factor  to  be  evaluated  over  a  wide  range 
of  crack  lengths.  Table  6  gives  the  values  of  the 
coefficients  of  the  polynomial  for  K  of  each  geometry. 

Figure  34  shows  the  variation  of  the  ratio  of  the 
stress  intensity  factor  to  nominal  plate  stress,  Op^,  as  t^ie 
crack  length  increased.  The  five  joint  conditions  shown  in 
Fig.  34  correspond  to  different  weld  sizes.  The  main  and 
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cross  plate  thickness  were  equal  and  the  weld  angle  is  equal 
to  45  degrees.  The  large  variation  of  the  stress  intensity 
factor  with  weld  si2e  is  readily  apparent.  A  decrease  in 
weld  size  resulted  in  a  substantial  increase  in  the  ratio 
of  K/ct  at  all  crack  sizes. 

P  Li 

Figure  35  shows  the  variation  in  stress  intensity 
with  weld  angle.  The  weld  cross  sectional  area  is  the  same 
for  each  weld  angle.  The  variation  in  the  stress  intensity 
with  weld  angle  is  seen  to  be  much  less  than  the  variation 
predicted  by  changing  the  weld  size  (H/Tp)  (see  Fig.  34). 
Increasing  the  weld  angle  from  30°  to  50°  results  in  approx¬ 
imately  a  ten  percent  increase  in  K/opi .  This  change  in 
weld  angle  represents  a  change  in  H./T  from  0.653  to  0.939  . 

A  similar  variation  in  II/Tp  in  Fig.  34  for  velds  with  equal 
legs  and  increasing  weld  size  would  decrease  the  stress 
intensity  factor  by  approximately  forty- five  percent.  The 
effect  of  weld  angle  on  the  stress  intensity  factor  for 
welds  of  equal  area  is  therefore'  small. 

The  similarity  in  the  shape  of  the  stress  intensity 
factor  curves  suggests  that  a  functional  relationship  may 

exist  which  describes  the  behavior  for  all  joint  qeomctries. 

(28) 

Harrison'  suggested  that  the  following  expression  could 
be  used  to  estimate  the  stress  intensity  factor  at  the  root 
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of  load  carrying  fillet  welds 

K  =  aw  /2w 1  tan  Tra/2w'  (21) 

aw  was  defined  as  the  nominal  gross  stress  acting  on  the  weld 
and  was  equal  to 

aw  =  P/B  (Tp  +  2H)  (22) 

This  estimate  considers  the  joint  to  be  equivalent  to  a 
plate  of  breadth  equal  to  Tp  +  2H  containing  a  central  crack 
of  length  2a. 

An  equation  similar  to  Eq.  22  was  found  to  provide 
a  reasonable  estimate  of  the  relationship  between  the  stress 
intensity  factor  and  crack  length.  The  more  accurate  secant 
finite  width  correction  was  used  in  place  of  the  tangent 
correction  used  by  Harrison.  In  addition,  a  non-dimensional 
correction  factor  was  developed  to  account  for  the  difference 
between  the  cruciform  joint  and  the  constant  breadth  plate. 
This  resulted  in  the  following  relationship  for  K  at  the 
weld  root: 

K  =  ct  A  /ifa  sec  it  a '  /2w '  (23) 

W  K 

where  A  is  non-dimensional  correction  factor  for  the 
K 

cruciform  shape. 

The  value  of  A  was  observed  to  be  a  function  of 
R 

H/Tp ,  the  weld  angle,  the  ratio  of  the  thickness  of  the  main 
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and  cross  plate,  and  the  value  of  a/w' .  The  value  of  X 

R 

was  calculated  for  each  geometry  by  equating  Eq.  23  to  the 

polynomial  expression  for  K  and  solving  for  X  .  The  heavy 

lines  in  Fig.  36  show  the  variation  in  X  with  a/w '  for 

R 

different  values  of  H/Tp. 

A  linear  approximation  was  fitted  to  the 
values  of  correction  factor  using  a  least  square  fit.  The 
fitted  equation  was 

X r  =  +  C2  a/w’  (24) 

The  linear  approximation  is  also  plotted  in  Fig.  36.  The 
linear  approximation  for  X  yielded  values  of  K  within  two 
percent  of  the  values  calculated  from  the  polynomial 
expressions  when  substituted  into  Eq .  23.  The  average  error 

was  less  than  one  percent  and  the  maximum  error  occurred  as 
a/w'  approached  :mro. 

The  coefficients  and  C9  were'  examined  to  find 
if  they  varied  systematically  with  joint  geometry.  Figure 
37  shows  the  variation  of  as  a  function  of  the  ratio 
H/T  .  It  is  apparent,  that  a  continuous  functional  relation¬ 
ship  exists  between  C|  and  II/Tp  over  a  wide  range  of  plate 
thicknesses  when  the  main  and  cross  plates  arc  of  equal 
thickness  and  weld  annie  is  4 X  degrees.  The  va  ues  of 
for  other  weld  ancles  are  also  shown  for  one  plate  thickness 
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and  a  constant  weld  area.  A  large  variation  in  with 
weld  angle  is  evident. 

Only  a  single  geometry  with  unequal  plate  thickness 
was  evaluated  since  only  one  series  of  tests  was  available 
in  the  literature.  This  condition  is  indicated  by  the  open 
circle  with  a  horizontal  line  in  Fig.  37.  The  value  of 
for  this  geometry  is  seen  to  be  slightly  less  than  a  joint 
with  equal  plate  thicknesses  for  the  same  value  of  H/Tp. 

Figure  38  shows  the  variation  of  C2  with  H/Tp.  It 
is  again  apparent  that  a  continuous  relationship  exists 
between  C2  and  H/Tp  for  a  weld  angle  of  45  degrees  and  Tp/Tc 
equal  to  unity.  As  observed  with  C^,  a  change  in  weld  angle 
results  in  a  substantial  change  in  C2 .  The  value  of  C2  for 
the  specimen  with  unequal  plate  thicknesses  is  seen  to  be 
greater  for  the  same  ratio  of  H/Tp. 

The  variation  of  and  with  weld  angle  as 

shown  in  Figs.  37  and  38  exaggerates  the  influence  of  the 

weld  anqle  on  the  stress  intensity  factor.  As  was  shown  in 

Fig.  35,  the  actual  variation  in  K/opL  with  weld  angle  was 

small.  The  large  variation  of  C-^  and  C2  is  due  to  the  fact 

that  K  is  expressed  as  a  function  of  <7  in  Eq.  23.  Since 

w 

the  weld  area  was  maintained  constant  in  the  analytical 
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study,  an  increase  in  the  weld  angle  decreases  the  nominal 
gross  weld  stress,  since  11  is  increased. 

Generally  the  load  carrying  capacity  of  a  joint 
is  expressed  in  terms  of  the  plate  force.  The  variation  of 
K  with  weld  angle  is  small  when  expressed  as  a  function  of 
the  plate  stress  as  was  shown  in  Fig.  35. 

The  values  of  and  C2  for  weld  angles  of  45 
degrees  and  Tp/Tc  equal  to  unity  were  fitted  by  least 
squares  to  polynomials  of  H/Tp.  This  yielded  the  following 
expressions : 

C,  =  0.52807  +  3.2872  (H/T  )  -  4.3610(H/T  )2 

^  ^  (25) 

+  3.6958  (H/Tn)  3  -  1.8745  (H/T  ) ^  +  0. 41495  (H/T_ ) 5 

r  r  r 

and 

=  0.21800  +  2.7173 (H/T  )  -  10.171(H/T  )2 

345  (26) 

+  13.  122  (Il/Tp)  J  -  7.754  6  <II/Tp)  H  +  1.7827 (H/Tp)  . 

These  equations  provide  estimates  of  and  C2  for  H/Tp 

between  0.20  and  1.20.  The  stress  intensity  factor  at  the 

weld  root  can  be  estimated  by  the  expression 

K  =  <J  [C.  +  C_(a/w')i  v  ■'ra  sec  tt1i,/2w'  (27) 

w  1  z 

for  joints  with  equal  plate  thicknesses  and  weld  anqles  of 
45  degrees. 
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The  value  of  K  predicted  by  Eq.  27  is  20  to  70 
percent  greater  than  the  estimate  provided  by  Harrison's 
approximation  (Eq.  21).  Equation  27  estimates  K  values 
within  two  percent  of  the  values  predicted  by  the  polynomial 
expressions.  The  value  of  K  for  joints  with  equal  plate 
thicknesses  and  other  weld  angles  (0  jt  45°)  can  also  be 
estimated  from  Eq.  27.  The  variation  of  K/qpi  for  various 
weld  angles  and  constant  weld  area  was  observed  to  be  small. 

A  reasonable  approximation  of  K  can  be  found  by  using  a  weld 
height  equivalent  to  a  45  degree  weld  of  the  same  area.  The 
equivalent  weld  height,  H,  can  be  found  from  the  relationship 
H  =  /H^W*  where  H*  and  W*  are  the  legs  of  the  actual  weld. 

Equation  27  can  be  expressed  in  a  more  useful 
design  form  as 

a 

K  =  (14.2h/t  )  tc1+C2  (a/w* )  J  /ira  sec  Tra/2w'  (28) 

where  ct  is  the  main  plate  stress.  This  equation  provides 
a  reasonable  estimate  of  the  stress  intensity  factor  at  the 
weld  root  in  terms  of  weld  and  plate  geometry  and  the  nominal 
plate  stress. 

4.3.2  Stress  Intensity  Factors  for  Toe  Cracks  of  Load 
Carrying  Fillet  Welds 

Table  7  summarizes  the  geometric  variables  for 
the  twenty- two  load  carrying  fillet  welded  joints  that  were 
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available  in  the  literature  and  this  study.  All  joints  had 
main  and  cross  plates  of  equal  thickness.  The  effect  of 
weld  angle  was  also  studied.  The  average  size  reported  in 
the  literature  for  each  series  was  used.  The  weld  size 
parameters  H  and  W  were  calculated  from  the  average  weld 
area  of  the  specimens.  A  comparable  method  was  used  for 
the  root  crack  analysis. 

The  influence  of  the  weld  penetration  was  also 
included  in  the  analysis.  The  average  depth  of  root 
penetration  for  each  geometry  is  given  in  Table  7.  The 
penetration  is  expressed  as  a  percentage  of  the  half 
thickness  of  the  main  plate.  Thirteen  geometric  conditions 
were  analyzed  with  variable  weld  penetration.  The  effect 
of  weld  penetration  was  included  in  the  analysis  since  a 
preliminary  study  had  indicated  that  the  stress 
at  the  weld  toe  was  a  function  of  the  weld  penetration. 
Increasing  the  weld  penetration  decreased  the  stress  at  the 
weld  toe.  Thirteen  of  the  twenty-two  geometries  were  an¬ 
alyzed  with  full  and  approximately  zero  weld  penetration  to 
determine  its  effect.  The  other  remaining  geometries  were 
analyzed  using  the  average  penetration  they  were  fabricated 
with . 


The  cracks  at  the  toe  of  the  weld  were  assumed  to 
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be  inclined  from  the  vertical  at  the  angle  of  13  degrees. 
This  angle  corresponded  to  the  condition  observed  during 
the  experimental  phase  of  this  investigation.  This  angle 
was  also  compatible  with  the  inclination  of  the  principle 
stress  axis  at  the  weld  toe. 

The  geometric  constant  w'  in  the  compliance 
polynomial  was  defined  as 

w'  =  Tp/  [ 2  cos  13°)  (29) 

The  compliance  of  each  joint  was  determined  for  thirteen 
crack  lengths  including  the  uncracked  joint.  The  ratios  of 
crack  length  to  w'  included  in  the  analysis  varied  from  0.15 
to  0.70  in  increments  of  0.05. 

A  sixth  order  compliance  polynomial  was  found  to 
provide  the  minimum  standard  error  of  estimate.  The 
standard  error  of  estimate  was  less  than  one-tenth  of  a 
percent  of  the  measured  compliances.  The  coefficients  of 
the  polynomial  expressions  for  the  stress  intensity  factor 
(see  Eq.  19)  for  each  geometry  is  given  in  Table  8. 

Figure  39  shows  the  predicted  variation  of  K/cFpj 
with  a/w'  for  two  different  plate  thicknesses.  The  effect 
of  weld  size  for  full  and  small  penetration  welds  is  also 
shown.  The  variation  of  H/Cp^  for  joints  with  full 
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penetration  welds  was  not  affected  by  weld  size.  However, 

for  joints  with  small  penetration  the  value  of  K/apL  is  seen 

to  increase  as  the  weld  size  decreased.  As  the  weld  size 

approached  the  plate  thickness,  the  effect  of  penetration 

on  the  ratio  K/a  decreased.  Partial  and  full  penetration 

PL 

welds  provided  about  the  same  conditions  at  the  weld  toe  for 
large  welds  (H/Tp  -  1). 

An  equation  comparable  to  Eq.  23  was  used  to 
express  the  stress  intensity  factor  in  terms  of  the  applied 
stress,  crack  size  and  plate  geometry 

K/ctPL  *  XT  /"3_  sec  ira/2w'  ( 30) 

where  w'  is  defined  in  Eq.  29  and  XT  is  a  non-dimensional 
correction  factor  that  is  analogous  to  the  stress  concen¬ 
tration.  The  form  of  Eq.  30  was  based  on  the  known  rela¬ 
tionship  for  a  finite  width  plate  with  a  through  crack.  This 
relationship  is  an  approximate  solution  for  a  double  edge 
crack  plate. 

Equation  30  was  equated  to  the  polynomials  that 
defined  K  as  a  function  of  a/w ’.  This  permitted  the 
variation  of  XT  with  crack  length  and  geometry  to  be 
determined.  Figure  40  shows  the  scatter  band  of  the 
variation  of  XT  with  a/w'  for  all  the  geometries  analyzed 
with  45  decree  weld  angles  and  full  root  penetration.  The 
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maximum  difference  between  XT  for  all  these  geometries  was 
less  than  51  for  all  crack  sizes.  The  magnitude  of  XT  was 
determined  at  values  of  a/w'  equal  to  0.02  through  0.10 
in  increments  of  0.02  and  from  0.10  through  0.70  in 
increments  of  0.05.  These  values  were  then  fitted  with  a 
polynomial  of  a/w'.  The  polynomial  which  provided  the 
minimum  standard  error  of  estimate  was 

XT  =  2.164  -  6.957  (a/w)  +  17.09  (a/w) 2  -  20.10  (a/w) 3 

+  8.842  (a/w)4  (31) 

This  equation  is  plotted  in  Fig.  40  as  a  solid  line  and  is 
seen  to  provide  a  reasonable  estimate  of  the  value  of  XT. 

Figure  41  shows  the  variation  of  XT  with  weld 
efngle.  The  weld  area  was  constant  for  each  joint  and  each 
joint  had  full  root  penetration.  Increasing  the  weld  angle 
is  seen  to  increase  XT  for  small  cracks.  This  effect 
diminishes  as  the  crack  grows.  As  a/w'  approaches  0.30  there 
is  no  significant  difference  due  to  the  angle  of  the  weld. 

Figure  42  shows  the  variation  of  XT  with  a/w'  for 
joints  with  small  root  penetration.  The  value  of  XT  is  seen 
to  increase  with  decreasing  weld  size.  This  behavior  was 
found  for  all  the  geometries  analyzed.  The  stress  intensity 
factor  at  the  weld  root  also  increases  with  smaller  weld 
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sizes  and  may  represent  a  more  critical  condition  when  the 
weld  penetration  is  small. 

A  stress  analysis  was  performed  on  selected  joints 
to  determine  if  the  behavior  of  the  stress  intensity  factor 
could  be  related  to  the  stress  concentration  at  the  weld  toe. 
Figure  43  shows  the  finite  element  mesh  used  in  the  stress 
analysis.  Constant  strain  triangular  elements  were  used. 

The  method  of  solution  is  given  in  Appendix  A.  The  stress 
at  the  weld  toe,  ^T0E/  was  taken  as  the  maximum  principle 
stress  in  the  shaded  element  (Fig.  43).  Figure  44  shows  the 
variation  of  the  stress  concentration  at  the  weld  toe 
apL)  with  H/Tp  for  full  and  zero  penetration  welds.  The 
plotted  points  are  from  joints  with  weld  angles  of  45  degrees. 
The  variation  of  the  stress  concentration  with  II/Tp  for 
joints  with  full  penetration  welds  is  small.  This  is  the 
same  behavior  that  was  found  for  XT  with  full  penetration 
joints.  Joints  with  zero  penetration  yielded  a  dramatic 
increase  in  the  stress  concentration  with  decreasing  weld 
size.  The  values  of  the  stress  concentration  for  zero 
penetration  welds  also  seem  to  lie  on  the  same  curve  for  all 
the  plate  thicknesses  included  in  the  analysis. 

The  variation  in  the  stress  concentration  at  the 
weld  toe  was  investigated  further.  Figure  45  shows  the 
variation  of  the  stress  concentration  factor  with  weld  angle 
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for  constant  area  welds.  Increasing  the  weld  angle 
increases  the  stress  concentration  for  joints  with  full 
and  zero  penetration  welds.  The  joints  with  full  penetration 
welds  have  lower  stress  concentrations.  The  variation  of 
the  stress  concentration  factor  with  the  weld  angle  is 
approximately  linear  for  both  weld  conditions. 

Figure  46  shows  the  variation  of  the  stress 
concentration  with  root  penetration  for  four  values  of  H/Tp. 
The  stress  concentration  for  small  welds  is  seen  to  increase 
much  more  rapidly  with  decreasing  weld  penetration  than  the 
larger  welds. 


The  value  of  the  stress  concentrations  at  the  weld 
toe  were  compared  to  the  values  of  XT  as  the  ratio  a/w' 
approached  zero.  This  limiting  value  of  XT,  which  was 
labeled  ~T,  was  found  by  taking  the  limit  of  Eq.  30  as  a/w' 
approached  zero.  This  yields 

K  =  apL  At  /va  (32) 


The  polynomial  expression  for  K  as  a/w'  approaches  zero 


yields 

T 

K  =  a  T  - [A,  (a/w' )  ]  2 

PL  2/wr  1 


Equating  Eqs.  32  and  3  3  provides  an  estimate  for  A^  as 


X 


(33) 


(A1/7r]iTp/[2w'] 


(34) 
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Figure  47  compares  XT  with  the  stress  concentration  at  the 
weld  toe  for  joints  of  various  plate  thickness,  weld 
penetrations,  and  weld  angle.  The  values  of  X"T  are 
compatible  with  the  stress  concentration  factor  although 
they  appear  to  be  slightly  larger.  Hence  the  stress  concen¬ 
tration  factor  provides  a  reasonable  estimate  of 

This  correlation  seems  reasonable  since  the  stress 

intensity  factor  for  small  cracks  at  the  weld  toe  should  be 

a  function  of  the  local  stress  field.  The  toe  stress 

from  the  finite  element  at  the  weld  toe  provides  an  estimate 

of  this  local  stress.  Since  a  rather  coarse  mesh  was  used 

to  estimate  the  stress  concentration,  (element  area  ^ 

2 

(w'/10)  )  the  true  stress  state  was  underestimated.  A  better 

estimate  of  the  stress  concentration  factor  would  provide 
better  agreement  with  X  .  The  dashed  line  in  Fig.  47 
represents  a  least  squares  fit  of  limiting  value  of  XT  and 
the  stress  concentration  factor  predicted  by  the  finite 
element  analysis.  \T  can  be  estimated  by 

XT  *  (aTOE/aPL)  +  °'2  (35) 

for  the  stress  concentrations  estimated  in  this  study.  If 
more  refined  methods  of  calculating  the  stress  are  used  the 
value  of  the  stress  concentration  should  approach  X^. 
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The  correlation  of  the  stress  concentration  factor 
with  ~T  provides  a  method  to  estimate  the  stress  intensity 
factor  K  for  various  crack  lengths  and  geometric  conditions 
which  were  not  considered  in  the  compliance  analysis. 

Figure  48  shows  the  variation  of  the  ratio  of  XT 
and  the  limiting  value  XT  with  a/w' .  The  relationships 
correspond  to  joints  with  45  degree  weld  angles  and  small 
weld  penetrations.  For  small  cracks,  a/w'  less  than  0.20, 
the  ratio  is  seen  to  be  approximately  the  same  for  all 
values  of  H/Tp.  The  curves  tend  to  diverge  at  larger  crack 
sizes. 


For  fatigue  crack  growth  the  estimate  of  K  for 
small  cracks  is  more  critical  since  the  largest  portion  of 
the  life  is  consumed  in  propagating  small  cracks.  An 
attempt  was  made  to  correlate  the  stress  intensity  factors 
at  the  weld  toe  for  both  full  and  partial  penetration  in 
terms  of  *T/*T  and  a/w '  . 

Equation  31  was  modified  by  expressing  the  ratio 
At/Xt  as  a  function  of  a/w'.  This  permitted  XT  to  be 
expressed  in  terms  of  its  value  as  a/w'  approaches  zero  (XT) 
for  various  crack  lengths.  Equation  31  in  this  form  yields 
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\T/XT  =  l-3.215(a/w')  +  7.897  (a/w 1 )  2  -  9.288(a/w')3 

+  4.086(a/w')4  (36) 

The  dashed  line  in  Fig.  48  compares  this  polynomial  with  the 
predicted  values  of  Equation  36  is  seen  to  provide 

a  reasonable  estimate  of  the  XT/X'T  ratio  for  these 
joints. 


The  stress  intensity  factor  for  cracks  at  the  toe 
of  load  carrying  fillet  welds  with  partial  root  penetration 
can  be  estimated  from: 

K  =  cxpL  XT  [l-3.215(a/w')+7.897  (a/w ' )  2-9.288(a/w' )  3 

+  4 . 086  (a/w' )  4]  /it a  sec  ira/2wl  (37) 

X^  can  be  estimated  from  the  stress  concentration  at  the 
weld  toe.  The  magnitude  of  the  stress  concentration  for  a 
particular  geometry  can  be  estimated  from  Figs.  45  and  46. 
Equation  37  is  accurate  within  21  for  full  penetration  welds 
with  45  degree  angle.  The  average  value  of  XT  for  all  the 
joints  analyzed  with  full  penetration  45  degree  welds  was 
found  to  be  2.16.  This  value  is  a  reasonable  estimate  of 
all  joints  with  full  penetration  welds.  The  accuracy  for 
joints  with  smaller  weld  penetration  is  less.  For  crack 
lengths  less  than  0.30w',  Eq .  37  provides  an  estimate  within 
5%  of  the  values  determined  from  the  polynomials. 
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4.3.3  Stress  Intensity  Factors  for  Toe  Cracks  of  Non- 
Load  Carrying  Fillet  Welds 

Test  data  were  available  for  thirteen  non-load 
carrying  weld  joints.  The  joint  geometries  are  summarized 
in  Table  9.  Specimens  NT-1  through  NT-9  were  evaluated 
based  on  the  average  dimensions  of  the  fatigue  specimens. 

The  weld  angle  was  varied  from  30  to  53  degrees  with  a 
constant  weld  area.  Specimen's  NT-10  geometry  was  based 
on  the  average  size  of  the  load  carrying  weld  fatigue 
specimen.  Specimens  NT-14  through  NT-19  were  selected  to 
provide  an  indication  of  variation  due  to  a  wide  range  of 
geometric  parameters. 

The  weld  root  penetration  was  taken  to  be  zero  for 
all  joints.  The  load  carrying  weld  joints  with  full  pene¬ 
tration  provide  information  that  would  be  directly  applicable 
to  the  non-load  carrying  joints  with  full  penetration. 

The  parameter  w'  in  the  compliance  polynomial  was 
defined  in  the  same  manner  for  load  carryinq  and  non-load 
carr,  inq  joint.'  ro  crack  growth  at  the  toe  (see  Fq .  20). 

This  assumed  the  crack  plane  to  be  13  degrees  from  the 
vertical.  The  same'  increments  of  a/w '  were  used  for  the 
compliance  ana  lj  ^  i  r.  of  both  weld  typos. 
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A  sixth  order  polynomial  was  again  found  to 
provide  the  minimum  standard  error  of  estimate.  The 
standard  error  of  estimate  was  less  than  one-tenth  of  a 
percent  of  the  measured  compliances.  The  coefficients  for 
the  stress  intensity  polynomials  for  each  joint  geometry 
are  summarized  in  Table  10. 

Figure  49  shows  the  typical  variation  of  K/apI 
with  a/w' .  The  value  of  K /anT  is  seen  to  increase  slightly 
with  an  increase  in  the  cross  plate  thickness.  Figure  50 
shows  the  variation  of  K/aOT  with  weld  angle.  The  scale 
of  the  abscissa  has  been  enlarged  to  show  the  variation  for 
small  values  of  a/w' .  The  value  of  K/apL  is  seen  to  in¬ 
crease  with  increasing  weld  angle.  The  sensitivity  of  the 
value  of  K  to  weld  angle  is  about  the  same  order  magnitude 
that  was  found  in  the  load  carrying  joints.  The  effect  of 
weld  angle  diminishes  with  increasing  crack  lengths.  This 
behavior  with  weld  angle  was  also  noted  in  the  load  carrying 
weld  joints. 

The  effect  of  weld  root  penetration  on  K/opL  is 
illustrated  in  Fig.  51.  The  full  penetration  joint  is  seen 
to  have  a  substantially  higher  value  of  K/op^  than  the 
penetration  joint  at  all  crack  sizes. 


zero 
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within  the  range  of  the  values  plotted.  The  value  of 
VXT  is  seen  to  be  approximately  the  same  for  all 
geometries  for  a/w'  less  than  0.25.  Some  divergence  is 
apparent  for  larger  cracks  for  the  different  ratios  of  H/Tp. 

Equation  36  which  was  developed  from  the  analysis 
of  load  carrying  joints  with  toe  cracks  is  plotted  in 
Fig.  54.  This  estimate  of  ^T/^T  is  seen  to  provide  a 
conservative  estimate  of  the  variation  of  *T/^”T  when  a/w' 
is  less  than  0.35.  The  value  for  larger  crack  sizes  falls 
within  the  band  of  the  values  for  the  various  ratios  of 
H/Tp.  Equation  36  provides  a  reasonable  estimate  of  the 
variation  of  ^T/^T  for  toe  cracks  in  both  load  carrying 
and  non-load  carrying  joints. 

The  value  of  XT  was  found  to  correlate  reasonably 
well  with  the  stress  concentration  factor  for  load  carrying 
joints.  No  such  correlation  was  found  for  the  non-load 
carrying  joints.  The  stress  concentration  factor  increased 
with  H/Td  for  joints  with  equal  plate  thicknesses,  T  /T  =1, 
and  with  the  weld  angle.  This  is  similar  to  the  behavior 
which  was  observed  for  load  carrying  joints  (see  Figs.  44 


and  45)  . 
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Attempts  to  correlate  the  stress  concentration 
factor  at  the  weld  toe  with  were  unsuccessful  as  is 
apparent  in  Fig.  55.  The  stress  concentration  is  seen  to 
decrease  with  increasing  XT.  The  stress  concentration 
overestimates  the  value  of  XT  for  most  of  the  geometric 
conditions  examined. 

The  stress  intensity  factor  for  toe  cracks  in 
non-load  carrying  joints  for  the  geometries  analyzed  can  be 
estimated  from  Eq.  37.  The  value  of  XT  for  these  geometries 
can  be  found  using  Eq.  34  with  the  leading  term  of  the 
stress  intensity  polynomial  given  in  Table  11.  The  value 
of  XT  for  other  similar  geometries  can  be  estimated  from 
these  values.  When  estimating  values  of  XT,  consideration 
should  be  given  to  the  fact  that  XT  increased  with  smaller 
welds  and  increasing  weld  penetration  increased  the  value 
of  ~T 
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5.  COMPARISON  OF  TEST  RESULTS  WITH  BEHAVIOR  PREDICTED 
BY  THE  FRACTURE  MECHANICS  OF  STABLE  CRACK  GROWTH 


The  fatigue  behavior  of  the  joints  tested  in  this 
study  as  well  as  the  results  of  previous  investigations  were 
evaluated  using  the  principles  of  the  fracture  mechanics  of 
stable  crack  growth.  The  equations  developed  from  the  simple 
power  law  of  crack  growth  were  used  to  relate  the  fatigue 
life  of  the  test  specimens  to  the  applied  stress.  The 
influence  of  joint  geometry  and  location  of  the  failure  were 
included  in  the  analysis  using  the  relationships  for  the 
stress  intensity  factor  developed  in  this  dissertation. 

5. 1  Crack  Growth  Behavior  of  Structural  Steels 

Figure  56  summarizes  the  results  from  five  crack 
growth  investigations  and  defines  the  crack  growth  rate, 
da/dN,  as  a  function  of  the  range  of  the  stress  intensity 
factor,  AK.  ^ )  The  growth  rates  are  for 

materials  with  yield  strengths  between  24  ksi  and  125  ksi. 
Also  included  are  data  from  Ref.  49  where  the  growth  rate  was 
measured  in  the  weld  metal  and  heat  affected  zone.  The  weld 
metal  yield  strengths  varied  from  70  to  90  ksi.  The  welds 
were  made  usinq  the  manual  shielded  metal  arc  process  and 
the  gas  metal  arc  process. 
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The  data  is  seen  to  fall  in  a  narrow  scatterband. 

The  slope  of  the  scatterband  is  about  one-third.  This  slope 

corresponds  to  an  exponent,  "n",  equal  to  3  in  Eq.  2. 

Barsom  and  Klingerman  observed  that  an  exponent  of  three 

provided  the  best  fit  to  their  test  data  when  the  power  law 

was  used  to  describe  the  relationship  between  growth  rate 

and  the  range  of  the  stress  intensity  factor.  Gurney  and 

Maddox  found  that  the  exponent  varied  slightly  with  material 

strength.  Their  studies  showed  a  variation  in  "n"  between 

2.4  and  3.4.  A  fatigue  study  of  the  effect  of  weldments 

on  steel  beams  has  shown  that  the  relationship  between  the 

stress  range  at  the  detail  and  its  life  yielded  a  value  of 

"m"  close  to  three  for  the  exponential  stress-life  model,  in 
(25) 

Eq.  8.  The  equivalence  of  the  exponent  "n"  in  the 

power  law  of  crack  propagation  to  the  exponent  "m"  in  the 
exponential  stress  life  relationship  was  discussed  in  part 
three.  This  assumption  was  tested  by  Hirt  in  his  analysis 
of  the  behavior  of  welded  beams.  He  found  that  an 

exponent  of  three  provided  a  reasonable  estimate  of  both 
"n"  and  "m"  in  Eqr, .  2  and  8  respectively.  An  exponent 
equal  to  1  was  also  selected  for  this  fatigue  cudy. 

Little  information  is  available  about  the  crack 
growth  rate  at  the  extreme  values  of  AK.  At  low  levels  of 
AK  a  threshold  in  crack  growth  rate  has  been  suggested.  This 
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has  been  investigated  by  Paris  in  ASTM  9310  steel. 

He  found  that  the  growth  rate  decreased  significantly  at 

AK  =  5.2  ksi  /Tn  .  Harrison  has  estimated  the  value  of  this 

threshold  level  from  an  analysis  of  rotating  beam  fatigue 
(52) 

results.  He  estimates  a  threshold  value  of  AK  between 

3.5  and  5.5  ksi  /Tn  for  steels.  These  estimates  are  shown 
in  Fig.  56  as  horizontal  lines  near  the  ordinate. 

At  higher  levels  of  AK,  Parry,  Klingerman  and 

Barsom  observed  that  the  crack  growth  rate  increased 
(37  46  47) 

markedly.  '  '  Parry  and  Klingerman  observed  this 

change  to  occur  when  the  net  section  stress  approached  the 
yield  point.  Barsom  suggested  that  the  increase  occurred 
at  a  critical  value  of  A6,  the  range  of  crack  opening- 
displacement.  The  bands  plotted  in  Fig.  56  also  include 
the  data  from  these  studies  at  high  AK  values.  The  rate  of 
crack  growth  observed  at  high  AK  levels  falls  within  the 
same  general  scatterband  as  the  lower  levels  of  AK.  The 
assumption  that  the  AK-da/dN  is  governed  by  Eq.  2  at  all 
levels  of  AK  seems  to  be  reasonable  in  view  of  the 
available  experimental  evidence. 

5 . 2  Analysis  of  Load  Carrying  Joints  Failing  from  the 
Weld  Root 


Six  different  joint  gecmetries  were  included  in 
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the  analysis.  Table  11  gives  the  average  dimensions  of  the 
joints  analyzed  and  their  references.  In  addition,  the 
analysis  included  data  from  a  crack  growth  experiment  on  a 
center  crack  plates  by  Klingerman.  The  dimensions  of 

the  plate  are  also  given  in  Table  11.  Tp  was  taken  as  the 
plate  width  and  the  initial  crack  length  is  expressed  as  a 
percentage  of  the  half  width  under  the  column  giving  the 
value  of  the  percentage  of  root  penetration.  This  data 
was  included  in  the  analysis  since  Eq.  27,  which  expresses 
the  value  of  K  for  root  cracks,  is  based  on  the  known 
relationship  for  K  for  a  center  cracked  plate.  By  equating 
the  coefficients  C2  and  the  weld  height,  H,  equal  to  zero 
and  C1=l,  Eq.  27  yields  the  K  value  for  the  center  cracked 
plate . 


The  plate  geometry  provides  a  means  of  evaluating 
the  error  in  the  estimate  of  the  stress  intensity  factor 
for  the  welded  joints.  The  stress  conditions  and  the 
geometry  of  the  plate  specimens  were  carefully  controlled 
thereby  minimizinn  any  associated  error  in  their  analysis. 

Equation  6  can  be  expressed  in  the  following  form 
for  an  exponent  of  3: 

-  (SR/3/F~  )"3  1/C 


N 


(38) 
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where  S_  refers  to  the  gross  section  stress  at  the  weld. 
R 

The  integral  "X"  is  defined  in  Eq.  5  as 


X  =  J  g  (a ,  w ' )  nda 

a  • 

1 


(39) 
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=  /[C1+C2  (a/w* )  ] /na  sec  7ra/2w '}  da 


a . 
1 


where  w'  is  defined  by  Eq.  20.  The  coefficients  and  C ^ 
are  defined  by  Eqs.  25  and  26  respectively  as  functions  of 
H/Tp. 


It  is  convenient  to  express  Eq.  38  in  terms  of  a 
"modified  stress  range"  defined  as 

SR*  =  SR/x1/3  (40) 

where  the  units  of  SR*  are  KSI  in  ^ ,  Substituting  Eq.  40 
into  Eq.  38  yields: 

N  =  Sr*"3/C  (41) 

This  equation  reduces  to  a  linear  form  with  the  following 
logaithnic  transformation 

log  N  =  -log  C  -  3  log  SR*  (42) 

The  load  carrying  joints  were  evaluated  by 
determining  S  *  for  each  joint.  The  integration  to 
determine  "X"  was  performed  numerically  using  twelve  point 
Gaussian  Quaduature.  Tests  of  the  integration  routine  for 
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comparable  integrals  have  yielded  results  to  within  0.1% 
of  the  closed  form  solution.  The  error  in  "X"  is  believed 
to  be  less  than  1%.  The  average  values  of  X^3  for  the  test 
specimens  analyzed  is  given  in  Table  11.  The  final  crack 
size,  a^,  was  taken  as  the  value  when  the  net  section  stress 
at  the  weld  was  equal  to  the  yield  stress.  The  influence  of 
the  final  crack  size  on  the  value  of  "X"  using  other 
approximations  such  as  the  ultimate  strength  caused  less  than 
a  2%  variation  in  "X".  The  crack  size  calculated  from  the 
yield  stress  criterion  gave  closer  agreement  to  the  measured 
size  in  the  fatigue  specimens  tested  in  this  study.  The 
value  of  a^  was  taken  as  half  the  thickness  of  the  main  plate 
minus  the  weld  penetration. 

Figure  57  summarizes  the  results  of  the  analysis. 
The  "modified  stress  range”  Sr*  is  plotted  as  a  function  of 
cycle  fife.  The  results  are  seen  to  fall  within  a  narrow 
scatter  band.  The  two  parallel  lines  shown  in  Fig.  57 
represent  the  variation  in  the  coefficient  "C"  found  by 
Barsom  in  his  study  of  the  crack  growth  behavior  of  ferrite- 
pearlite  steels.  These  lines  are  plotted  at  a  slope  of  -1/3, 
correspondinq  to  t  lie  exponent  n=3  in  the  crack  growth  power 
law. 


Tlie  results  from  K 1  ingerman '  s  cracked  plate  study, 


shown  as  the  open  squares,  are  seen  to  fall  within  the 


-66 


limits  for  the  welded  joints.  The  estimates  of  the  K 
relationship  for  welded  joints  developed  in  this  study  seem 
to  be  reasonable. 

The  large  scatter  in  the  data  reported  in  Ref.  15 
is  believed  to  be  due  to  the  use  of  average  values  of  the 
initial  crack  size,  a^,  and  weld  size.  The  welding  process 
was  the  same  for  all  the  specimens  plotted  from  Ref.  15. 
Since  wide  variations  in  the  average  values  existed  for  each 
joint  of  this  study,  comparable  variation  in  the  four  welds 
could  be  expected  for  the  joints  reported  in  Ref.  15.  The 
fatigue  failure  should  occur  at  the  weld  with  the  smallest 
penetration  (largest  a^)  and  size.  The  use  of  average 
geometry  and  crack  size  would  obviously  lead  to  increased 
scatter  as  observed  in  Fig.  57.  The  degree  of  weld 
penetration  and  consequently  the  initial  crack  size  was 
more  precisely  defined  by  the  other  investigators. 

The  linear  form  of  the  log  S  *  and  log  life 
relationship  is  seen  to  be  valid  over  a  wide  interval  of 
life.  Lives  from  10^  to  5x10*^  cycles  arc  seen  to  follow  the 
linear  relationship  predicted  by  Eg.  42.  The  equation  is 
seen  to  be  valid  in  both  the  high  stress-low  cycle  life  and 
low  stress-high  cycle  life  regions.  The  results  from  a 
fatigue  study  in  either  of  these  regions  can  be 
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extrapolated  over  a  wide  range  of  stress  and  fatigue  life. 

At  low  levels  of  S  *  and  high  cycle  life  the 
scatter  is  seen  to  increase.  The  estimate  of  AK  for  the 
average  initial  crack  size  for  joints  from  Ref.  15  tested 
at  SR*=11  is  7.0  ksi  /In  .  The  crack  growth  rate  at  this 
level  of  AK  has  not  been  established  experimentally. 

Figure  56  indicates  that  at  these  low  levels  of  AK,  the 
threshold  of  crack  growth  predicted  by  Paris  and  Harrison 
is  approached. 

The  fatigue  life  of  other  load  carrying  joints 
that  exhibit  failure  from  the  weld  root  can  be  estimated 
from  Fig.  57.  The  magnitude  of  "X"  can  be  estimated  by 
integrating  Eq.  39  between  the  limits  provided  by  the 
initial  crack  size,  a^,  and  the  final  crack  size,  a^..  The 
final  crack  size  can  be  estimated  as  the  crack  size  at 
which  the  net  section  stress  of  the  weld  is  equal  to  the 
yield  stress.  The  initial  crack  size,  a^,  is  a  function  of 
the  weld  root  penetration.  Therefore,  for  a  given  weld  size, 
plate  size  and  alternating  plate  load,  the  degree  of  weld 
penetration  necessary  to  provide  the  desired  fatigue  life 


can  be  determined. 
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5 . 3  Analysis  of  Load  and  Non-Load  Carrying  Weld 
Toe  Failures 


Three  non-load  and  three  load  carrying  fillet 
welded  joints  were  analyzed  for  crack  growth  at  the  weld 
toe.  The  joint  geometries  are  summarized  in  Table  12.  The 
analysis  of  crack  growth  at  the  fillet  weld  toe  was  analo¬ 
gous  to  the  study  of  root  cracks.  The  stress  intensity 
factor,  K,  was  estimated  from  Eq .  37.  The  equation  can  be 
expressed  in  the  following  form 

K  =  apL  XT  h (a , a/w 1 )  (43) 

where 

h(a,a/w')=  /ra  sec  iTa/2w'  (44) 

The  function  AT/~T  was  defined  by  Eq .  36  in  terms  of  a/w'. 
The  expression  for  the  fatigue  life  can  be  obtained  from 
Eq.  5  and  is 

N  =  (Tt  Sr)"3  Y/C  (45) 


The  initial  crack  size,  a.,  existing  at  the  too 

of  fillet  welds  has  been  examined  at  the  Welding 
(21  27) 

Institute.  '  The  toes  ol  fillet  welds  made  be  several 

welding  processes  were  examined  metalographically  by  taking 
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slices  through  the  thickness  at  the  toe.  These  studies 

confirmed  that  crack  like  instrusions  of  slag  and  other 

inpurities  existed  at  the  weld  toe.  These  intrusions  as 

well  as  crack  like  cold  laps  and  undercut  were  observed  to 

provide  initiation  sites  for  fatigue  crack  growth.  The 

depth  of  these  defects  were  found  to  range  from  0.0004  in. to 

0.02  in.  Welded  specimens  tested  with  artificial  notches 

cut  at  the  weld  toe  to  a  depth  of  0.005  in.  yielded  fatigue 

lives  comparable  to  those  observed  from  specimens  without 
(27 ) 

notches.  The  fatigue  life  was  shown  to  increase  when 

the  toes  were  treated  by  grinding  or  peening  to  remove  or 
reduce  the  initial  flaw.  Also  the  life  increased  when  the 
welding  process  was  modified  to  decrease  the  initial  flaw 
size.  The  fatigue  behavior  of  a  fillet  weld  joint  is 
consequently  a  function  of  crack  growth  since  the  initial 
crack  like  defects  exist  at  the  weld  toe. 

The  integral  Y  which  expresses  the  relationship 
between  the  fatigue  life  and  applied  stress  range  is  a 
function  of  a.,  a.  ami  w'.  The  value  of  Y  was  determined 
numerically  for  various  values  of  a^  ,  a^  and  w'.  Integration 
was  again  performed  by  the  twelve  point  Gaussian  Quadrature 
A  singularity  in  the  integrand  exists  as  a^  approaches  z  ro. 
The  integration  was  performed  by  dividing  the  crack  interval 
between  a.  and  af  into  four  parts.  The  first  increment  was 


taken  as  a^  to  2a^  to  account  for  the  singularity. 

Figure  58  shows  the  variation  of  Y  with  for 
two  plate  thickness.  The  final  crack  size,  a^,  was  also 
varied.  a ^  was  taken  as  0. 50w'  =  0.0641  in.  for  the  4-  in. 
plate  and  0.80w'=0.821  in.  for  2  in.  thick  plate.  These 
values  provide  an  extreme  variation  in  plate  thickness  anci 
final  crack  size.  Values  of  Y  for  other  plate  thickness  and 
final  crack  sizes  were  bounded  by  the  curves  shown  in 
Fig.  58.  The  variation  of  Y  with  plate  thickness  and  final 
crack  si2e  is  seen  to  bo  negligible  compared  to  the 
influence  of  the  initial  crack  size.  This  is  particularly 
true  for  small  crack  sizes.  At  larger  crack  sizes  the 
value  of  Y  was  approximately  constant.  Hence,  Y  is 
primarily  a  function  of  a^,  the  initial  crack  size. 

The  value  of  a ^  and  the  corresponding  value  of  Y 

was  not  estimated  directly  in  the  analysis.  The  fatigue 

behavior  of  joints  failing  from  t: h« ■  wold  too  was  estimated 

by  assuming  that  Y  was  the  sane  for  all  specimens.  Hence, 

the  cyclic  life  of  each  specimen  was  plotted  is  a  function 

of  the  parameter  X  S  .  Figure  r»‘>  summarizes  the  available 

i  K 

test  data.  The  data  shown  includes  both  load  and  non-load 
carrying  fillet  welded  specimens.  The  value  of  Y  for  each 
geometry  was  estimated  as  outlined  in  part  four.  The  value 
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of  XT  estimated  for  each  geometry  is  given  in  Table  12. 

The  scatter  bands  provided  by  the  test  data  are 
plotted  at  a  slope  of  -1/3.  The  data  all  fall  within  a 
reasonably  narrow  band.  The  use  of  a  slope  equal  to  -1/3 
is  seen  to  provide  a  reasonable  estimate  of  the  linear 
relationship  between  ~T  SR  and  the  cycle  life.  The  variation 
in  the  test  data  is  greater  than  predicted  for  the  root 
cracks  shown  in  Fig.  57.  This  increase  in  the  scatter  can  be 
attributed  to  the  variation  in  initial  flaw  size,  a^,  and  the 
variation  of  the  crack  growth  constant  for  the  steels  used  in 
the  fabrication  of  the  specimens.  Obviously  differences  in 
the  welding  processes  could  substantially  affect  the  initial 
crack  size,  a^.  The  quality  of  the  weld  for  a  particular 
welding  process  can  also  vary  due  to  uncontrolled  variables 
such  as  operator  performance.  Manual  shielded  arc  welds  were 
used  in  Refs.  11  and  18.  The  manual  welded  specimens  from 
Ref.  11  arc  plotted  as  open  circles.  These  data  define  the 
lower  limit  of  the  scatter.  The  upper  limit  of  the  data  is 
provided  by  manually  welded  specimens  from  Ref.  18.  These 
data  are  shown  by  the  closed  circles.  The  automatic  welds 
used  in  the  fabrication  of  the  specimens  of  the  experimental 
study  fall  near  the  center  of  the  scatter  bands. 


The  initial  crack  size  can  be  estimated  by  eval- 
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uating  the  ratio  Y/C  in  terms  of  the  scatter  bands  shown  in 

Fig.  59.  The  value  of  Y/C  corresponding  to  the  upper  band 

is  1.1x10^  and  the  lower  band  is  1.2x10^.  Barsom's  crack 

growth  data  provided  an  estimate  of  the  magnitude  of  C.  He 

observed  C  to  vary  between  1.37  to  3.6x10  for  ferrite- 

pearlite  steels.  The  corresponding  values  of  Y  vary 

between  1.65  and  39.5.  If  average  values  are  used,  Y/C  = 

6.1xl010,  C=2 . 5xl0~10 ,  then  Y  -  15.  These  values  of  Y  are 

compatible  with  the  measured  initial  crack  sizes  found  at 

(21  27) 

fillet  weld  toes  as  shown  in  Fig.  58.  A  value  of  15 

for  Y  corresponds  to  an  initial  crack  size  of  about  0.001  in. 
in  Fig.  58. 


The  scatter  for  the  results  from  an  individual 
study  are  seen  in  Fig.  59  to  be  much  less  than  the  total 
scatter  for  all  results  plotted.  The  dashed  lines  in  Fig. 

59  show  the  scatter  in  the  results  from  the  fatigue  tests 
from  this  study.  Consequently  the  values  of  Y  calculated 
from  the  overall  scatterbands  and  the  limits  of  the  C  values 
from  Barsom's  study  overestimate  the  variation  of  the  initial 
flaw  size  for  a  particular  test  series.  However,  unless  the 
initial  flaw  size  is  measured  for  each  specimen  or  an 
estimate  of  the  value  of  C  corresponding  to  the  material  used 
to  fabricate  the  specimen,  the  two  parameters  Y  and  C  are 

confounded.  An  estimate  of  the  variation  in  a.  for  the 

i 

results  of  the  present  study's  fatigue  specimens  was  made  by 
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taking  C  equal  to  3.6x10  This  value  of  C  was  estimated 

from  Barsom's  study  and  was  felt  to  be  a  reasonable  estimate 
for  A514  steel.  Using  this  value  for  C  and  the  values  of 
XT  Sr  and  cycle  life  defined  by  the  dashed  lines  in  Fig.  59 
yielded  a  variation  of  a^  between  0.0015  in.  and  0.007  in. 
These  values  of  a^  correspond  to  values  of  Y  of  12.6  and  7.31 
respectively.  They  provide  a  more  concise  measure  of  the 
variability  of  initial  flaw  size  at  the  weld  toe.  These  flaw 
sizes  corrode  ^o  equivalent  lengths  of  sharp  cracks  at  the 
weld  toe.  A  blunt  intrusion  at  the  weld  toe  of  greater  depth 
may  provide  a  longer  life  then  inferred  from  this  analysis 
since  fatigue  crack  initiation  may  consume  an  appreciable 
portion  of  the  life. 

The  linear  trend  between  \  SD  and  cycle  life  in 

Fig.  59  is  seen  to  be  valid  over  a  large  range  of  cycle 

lives.  The  range  of  lives  included  in  the  analysis  were 
3  6 

2x10  to  2x10  cycles  of  load.  These  lives  cover  the  full 
range  of  the  high  cycle  fatigue  and  into  the  low  cycle 
regime  of  fatigue. 

The  fatigue  life  of  a  load  and  non-load  carrying 
weld  joint  can  be  estimated  directly  from  Fig.  59  when  fail¬ 
ure  occurs  from  crack  growth  at  the  weld  toe.  The  factor  X^ 

can  be  estimated  from  the  relationships  developed  in  this 
stuck/ . 


r 
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6.  CONCLUSIONS 


The  main  findings  of  this  study  are  summarized  in 
this  chapter.  They  are  based  on  the  analysis  of  the  fatigue 
performance  of  fillet  welded  joints  stressed  perpendicular 
to  the  weld  line. 

1.  The  stress  intensity  factor  for  root  cracks  in 
load  carrying  joints  was  found  to  be  a  function 
of  weld  size,  weld  angle,  and  the  thickness  of 
the  cross  plates.  A  relationship  was  developed 
which  permitted  the  stress  intensity  factor,  K, 
to  be  estimated  for  equal  leg  welds  and  equal 
thickness  plates. 

2.  The  stress  intensity  factor  for  toe  cracks  in 
load  carrying  joints  was  found  to  be  proportional 
to  the  local  stress  at  the  weld  toe.  This  stress 
was  estimated  from  a  finite  element  analysis  and 
found  to  be  a  function  of  weld  angle,  weld  size, 
plate  thickness  and  weld  penetration.  Increasinq 
the  weld  penetration  decreased  the  stress  and  the 
K  value  as  did  increasinq  weld  size.  The  value 
of  K  for  full  penetration  welds  was  found  to  be 
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approximately  the  same  for  all  geometries.  The 
stress  intensity  relationship  was  developed  in 
terms  of  a  non-dimensional  correction  factor  AT> 
The  value  of  XT  was  observed  to  decay  as  crack 
length  increased.  XT,  the  value  of  XT  as  the 
crack  length  approached  zero,  was  found  to  be 
proportional  to  the  stress  concentration  at  the 
weld  toe.  The  decay  function  of  XT  as  crack 
length  increased  was  found  to  be  the  same  for  all 
geometries . 

3.  The  stress  intensity  factor  for  several  non-load 
carrying  joints  was  obtained.  The  relationship 
developed  to  predict  K  for  load  carrying  weld  with 
toe  cracks  was  found  to  be  applicable.  However 

no  correlation  was  found  between  the  stress  at 
the  weld  toe  and  the  stress  intensity  factor. 
Further  study  is  desirable.  The  value  of  K  at  the 
weld  too  was  observed  to  increase  as  the  weld 
penetration  increased  and  weld  size  was  decreased. 

4.  The  failure  mode  of  load  carrying  welds  at  high 
stress  levels  is  dependent  upon  the  weld  size  and 
the  heat  affected  zone  inclusion  content.  A 
shearing  mode  of  failure  along  the  fusion  lino  w^s 
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exhibited  by  joints  with  small  welds.  This 
behavior  was  attributed  to  cracking  along 
included  particles  at  high  stress  levels. 

5.  The  fatigue  behavior  of  both  the  non-load 
carrying  and  load  carrying  joints  was  analyzed 
using  fracture  mechanics  principles.  The 
fatigue  behavior  of  many  different  specimen 
geometries  were  correlated.  The  results  provide 

a  rational  means  of  estimating  the  fatigue  perfor¬ 
mance  of  a  joint  for  lives  between  10 3  and 
5xl06  cycles  of  load  application.  The  correlation 
provided  at  low  cycle  lives  provides  a  means  to 
determine  the  fatigue  performance  of  structural 
joints  subjected  to  high  stress  ranges  from  a 
knowledge  of  their  low  stress  performance. 

6.  The  influence  of  the  initial  crack  size  on  the 
life  of  toe  crack  failures  was  found  to  be  great. 
The  life  increased  rapidly  as  the 

initial  crack  sice  approaches  zero.  At  crack 
sizes  greater  than  0.01,  the  effect  was  much  less. 
Treatment  of  the  weld  toe  by  grinding  and  peening 
to  reduce  the  initial  crack  size  offer  means 


A 


r 
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of  substantially  increasing  the  fatigue  life. 

7.  The  threshold  value  of  AK  for  toe  cracks  appeared 
to  be  about  2  ksi  /Tn  ,  much  lower  than  estimated 
by  others.  The  value  for  root  cracks  appeared  to 
be  about  7  ksi  /Tn  .  The  fatigue  behavior  of 
welded  joints  at  low  stress  ranges  is  dependent 
upon  the  crack  growth  behavior  at  these  low  levels 
of  AK.  Crack  growth  studies  at  low  levels  of  AK 
are  needed  to  allow  the  accurate  prediction  of  the 
fatigue  performance  of  joints  at  low  levels  of 
stress . 

8.  The  finite  element  technique  of  stress  analysis 
was  shown  to  provide  a  method  of  estimating  the 
stress  intensity  factor  for  cracks  in  fillet 
welded  joints.  This  technique  provides  a 
powerful  and  straight  forward  means  to  estimate 
the  stress  intensity  factor  for  other  complex 
geometries.  The  analysis  of  the 

cracking  behavior  of  structures  by  fracture 
mechanics  need  no  longer  be  restricted  to  simple 
qeometrics  with  known  stress  intensity  factors 
since  these  may  be  estimated  by  a  finite  element 


analysis . 


TABLES 
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TABLE  1 

MECHANICAL  AND  PHYSICAL  PROPERTIES  OF  PLATF  MATFRTAL 
ASTM  GRADE  A514F 

COMPOSITION 

(PERCENTAGE) 

C  Mn  P  S  Cu  Si  Ni  Mo  _ V_  B 

0.15  0.85  0.010  0.018  0.17  0.17  0.75  0.45  0.06  0.0045 

MECHANICAL  PROPERTIES 
YIELD  STRENGTH  =108  KSI 

TENSILE  STRENGTH  =118  KSI 

ELONGATION  IN  2  IN.=  36% 


REDUCTION  IN  AREA 


54% 
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TABLE  2 


EXPERIMENT  DESIGN  FOR  FATIGUE  TESTS 


30 


NOMINAL  STRESS  RANGE  MAIN  PLATE 
KSI 


40 


50 


60 


70 


80 


90 


to 
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_1  C/3 
<  * 
55 

M  W 

T.  H 

S3 

O. 

|  55 

S  *2 


50 

60 

70 

80 

90 

100 

110 


3* 


♦Number  of  Specimens  Tested. 
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TABLE  3 

FATIGUE  TEST  RESULTS 
LOAD  CARRYING  WELD  JOINTS 


Specimen 

Stress  Range 
KSI 

Maximum  Stress 
KSI 

Fatigue  Life 
Cycles 

Failure 

L-34 

90 

110 

610 

Shear 

L-28 

90 

110 

1,548 

Root 

L-7 

90 

110 

1,515 

Shear 

L-19 

80 

100 

4,224 

Toe 

L-17 

80 

100 

4,726 

Toe 

L-39 

80 

100 

740 

Root* 

L-33 

80 

100 

2,425 

Toe 

L-l 

70 

90 

1,267 

Root* 

L-29 

70 

90 

3,147 

Root 

L-36 

70 

90 

3,737 

Shear 

L-5 

70 

90 

2,900 

Toe 

L-18 

60 

80 

6,948 

Root 

L-10 

60 

80 

7,497 

Root 

L-14 

60 

80 

6,681 

Root* 

L-6 

50 

70 

7,721 

Root 

L-26 

50 

70 

13,060 

Root* 

L-21 

50 

70 

17,839 

Toe 

L-4 

40 

60 

22,026 

Root 

L-9 

40 

60 

19,878 

Root 

L-20 

40 

60 

23,422 

Root 

L-23 

30 

50 

36,091 

Root 

L-24 

30 

50 

35,526 

Root 

L-ll 

30 

50 

42,959 

Root 

*Porous  Welds. 


TABLE  4 

FATIGUE  TEST  RESULTS 
NON-LOAD  CARRY  WELD  JOINTS 
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Specimen 

Stress  Range 
KSI 

Maximum  Stress 
KSI 

Fatigue  Life 
Cycles 

Failui 

N-6 

90 

110 

8,683 

Toe 

N-18 

90 

110 

6,852 

Toe 

N-35 

90 

110 

5,319 

Toe 

N-25 

80 

100 

16,550 

Toe 

N-7 

80 

100 

13,896 

Toe 

N-17 

80 

100 

13,370 

Toe 

N-10 

70 

90 

19,945 

Toe 

N-12 

70 

90 

22,107 

Toe 

N-33 

70 

90 

11,780 

Toe 

N-16 

60 

80 

42,331 

Toe 

N-37 

60 

80 

35,536 

Toe 

N-30 

60 

80 

22,205 

Toe 

N-8 

50 

70 

41,195 

Toe 

N-26 

50 

70 

61,858 

Toe 

N-4 

50 

70 

34,830 

Toe 

N-2I 

40 

60 

99,446 

Toe 

N-3 

40 

60 

58,321 

Toe 

N-9 

40 

60 

87,485 

Toe 

N-29 

30 

50 

225,878 

Toe 

N-ll 

30 

50 

217,998 

Toe 

N-39 

30 

50 

182,390 

Toe 

DIMENSIONS  OF  LOAD  CARRYING  JOINTS  ANALYZED  FOR  ROOT  CRACKS 
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STRESS  INTENSITY  FACTOR  POLYNOMIAL  COEFFICIENTS  FOR 


ROOT  CRACKS  OF  LOAD  CARRYING  WELD  JOINTS 

nA^ 


SPECIMEN 

2A1 

3A, 

4A3 

5A4 

LR  - 

1 

11.812 

24.199 

_ 

71.383 

121.45 

LR  - 

2 

11.671 

24.270 

- 

71.220 

121.21 

LR  - 

3 

10.856 

24.377 

- 

71.858 

122.38 

LR  - 

4 

10.930 

24.859 

- 

73.045 

123.26 

LR  - 

5 

10.690 

24.853 

- 

73.181 

123.53 

LR  - 

6 

15.495 

19.418 

- 

68.691 

121.62 

LR  - 

7 

15.535 

19.655 

- 

69.761 

122.67 

LR  - 

8 

17.655 

11.122 

- 

54.809 

114.20 

LR  - 

9 

9.1429 

26.497 

- 

81.572 

131.96 

LR  - 

10 

15.899 

18.088 

- 

66.556 

120.39 

LR  - 

11 

18.816 

4.8185 

- 

41.964 

106.31 

LR  - 

12 

6.8059 

32.931 

-106.76 

151.33 

LR  - 

13 

9.9418 

24.911 

- 

74.884 

125.69 

LR  - 

14 

12.078 

24.254 

- 

71.443 

121.15 

LR  - 

15 

13.831 

23.329 

- 

73.408 

123.55 

LR  - 

16 

15.345 

20.071 

- 

70.008 

122.50 

LR  - 

17 

16.631 

15.635 

- 

63.019 

118.91 

LR  - 

18 

8.5541 

27.607 

- 

86.785 

136.55 

LR  - 

19 

11.612 

24.112 

- 

70.835 

120.93 

LR  - 

20 

13.835 

23.200 

- 

72.980 

123.12 

LR  - 

21 

15.680 

19.120 

- 

68.723 

121.94 

LR  - 

22 
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TABLE  10 


STRESS  INTENSITY  FACTOR  POLYNOMIAL  COEFFICIENTS  FOR 
TOE  CRACKS  OF  NON-LOAD  CARRYING  JOINTS 


SPECIMEN 

2A1 

3A2 

nAj 

4A3 

5A4 

6A5 

NT 

1 

15.347 

-  73.968 

222.40 

-308.61 

170.38 

NT 

- 

2 

17.240 

-  89.506 

271.99 

-378.81 

207.07 

NT 

- 

3 

18.288 

-  98.451 

300.94 

-420.06 

228.70 

NT 

- 

4 

18.936 

-104.12 

319.47 

-446.59 

242.64 

NT 

- 

5 

19.824 

-112.12 

345.88 

-484.55 

262.64 

NT 

- 

6 

20.352 

-117.04 

362.29 

-508.25 

275.16 

NT 

- 

7 

20.827 

-121.58 

377.60 

-530.46 

286.90 

NT 

- 

8 

21.434 

-127.63 

398.26 

-560.59 

302.89 

NT 

- 

9 

21.912 

-132.69 

415.91 

-586.55 

316.72 

NT 

- 

10 

20.365 

-117.12 

362.53 

-508.59 

275.34 

NT 

- 

11 

14.475 

-  82.982 

275.24 

-360.02 

194.92 

NT 

- 

12 

17.086 

-  96.572 

295.32 

-402.93 

213.33 

NT 

- 

13 

14.964 

-  85.953 

265.97 

-372.09 

201.15 

NT 

- 

14 

25.678 

-144.86 

444.11 

-605.80 

321.46 

NT 

- 

15 

8.0191 

-  46.212 

142.83 

-200.35 

108.30 

NT 

- 

16 

11 . 300 

-  64.902 

200.80 

-281.14 

152.09 

NT 

- 

17 

45.063 

-275.68 

912.88 

-1318.9 

728.01 

NT 

- 

18 

14.730 

-  83.529 

254.31 

-346.75 

183.00 

NT 

- 

19 

13.361 

-  75.433 

230.87 

-314.88 

166.81 

AVERAGE  DIMENSIONS  OF  ROOT  CRACK  FATIGUE  SPECIMENS  ANALYZED 
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*Non  Load  Carrying  Welds-Zero  Penetration  Assumed 


(a)  Non- Load  Carrying 


(b)  Load  Carrying 


Nominal  Dimensions 

Tp  =  Tq  =  0.625  in. 

W  =  H  =0.500  in. 

Hc  =  2.00  in. 

L  =  24.0  in. 

Width  =  1.200  in. 

Fig.  1  Fatigue  Specimen  Dimension 
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Fig.  2  Welding  Sequence 
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Fig. ' 3  Etched  Cross  Section  of  Typical  Joint 
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Fig.  4  Toe  Failure-Non-Load  Carrying  Joint 


Fig.  5  Rogt  Failure-Load  Carrying  Joint 
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CO  JC 


CYCLES  TO  FAILURF 

Fig.  7  Fatigue  Life  of  Load  and  Non-Load  Carrying  Joints 
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Weld  Angle  -  Q j  =  Arctan  [^i/wj] 
Angle  of  Toe  Crack  =  a 


Fig.  8 


Definition  of  Weld  and  Crack  Angle 
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WELD  ANGLE 

Fig.  9  Histogram  of  Weld  Angle-Non-Load  Carrying  Joints 


9/ a 

°  FAILURE 


Fig.  10  Histogram  of  9  to  9pailurc 


-Non  Load  Carrying  Joints 
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Fig.  15  Toe  Crack  Growth  Along  Fusion  Line 


Fig.  16  End  of  Fusion  Line  Crack  -  66X 
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Pig.  19  Crack  Branching  to  Included  Particles  -  66X 
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Fig.  20  Crack  Growth  Path  at  Included  Particles  -  133X 
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Fig.  21  Root  Crack  Tip  Untested  Specimen  133X 
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Fig.  22  Root  Crack  Tip  Untested  Specimen  With  Larqo  Plate 
Gap  -  133X 
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Fig.  24  Transgranular  Root  Crack  Growth  -  133X 
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Fig.  25  Vortical  r’nd  of  Root  Crack  At  Failure  -  66X 


Fig.  26  Failure  of  Joint  With  i-orous  Weld 
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Load  Carrying  Joint  Root  Crack  Finite  Element  Mesh 


Load  Carrying  Joint  Toe  Crack  Finite  Element  Mesh 
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Fig.  32  Non- Load  Carrying  Toe  Crack  Finite  Element  Mesh 


Tp=  Tc  =  0.662  in. 
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Fig.  34 


Variation  of  K/a  .  With  !I/T  for  Root  Cracks  in 
Load  Carrying  Joints  ^ 
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Tp  in-  Tp/Jc 


Fig.  37  Variation  of  the  Coefficient  C.  With  H/T 
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of  the  Coefficient 


Fig . 
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40  Variation  of  for  Load  Carrying  Joints  with  Full 
Penetration  ana  45°  Degree  Welds 


r  Stress  Analysis  of  Load 


Penetration  |TP=TC 


1.260 


0.786 


0.764 


0.662 


0.630 


0.625 


Fig.  45  Change  in  Stress  Concentration  at  Weld  Toe  With 
Weld  Angle 


Fig.  46  Change  in  Stress  Concent  ration  at  Weld  Toe  With 
Root  Penetration 
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Fig.  47  Correlation  of  '  With  Stress  Cone 
Weld  Toe  ! 
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With  Crack  Length  -  Small 
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Fig.  53  Change  in  >  With  Weld  Penetration 
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APPENDIX 


Methods  of  Solution  of  Finite  Element  Equations 


The  method  of  solution  used  to  solve  the  set  of 
equations  resulting  from  a  finite  element  analysis  becomes 
important  when  the  size  of  the  system  of  equations  is  large. 
The  mesh  used  to  determine  the  stress  distribution  in  the 
joint  had  approximately  four-hundred  nodal  points.  This 
resulted  in  an  800  x  800  stiffness  matrix.  The  solution  for 
this  mesh  was  found  using  the  Gauss-Seidel  procedure  with 
over  relaxation.  The  form  of  the  program  and  its  logic  was 
based  on  the  program  in  Ref.  53.  This  method  was  used  since 
only  the  non-zero  terms  of  the  stiffness  matrix  need  to  be 
stored.  The  displacement  of  all  nodal  points  are  solved  and 
can  be  used  to  determine  element  stresses.  The  solution  was 
assumed  to  have  converged  when  the  percent  change  in  the 
displacement  at  the  loaded  nodes  at  the  end  of  the  main 
plate  was  less  than  one-tenth  of  a  percent. 


The  mesh  used  to  determine  the  compliance  of  the 
joints  was  considerably  more  refined.  Fifteen-hundred  nodal 
points  were  used  and  resulted  in  a  global  stiffness  matrix 
of  3000  x  3000.  The  finite  element  mesh  away  from  the  crack 
itself  does  not  change  with  crack  length  and  therefore  the 
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stiffness  matrices  of  these  areas  are  constant.  The  method 
of  solution  made  use  of  the  constant  stiffness  of  these  areas 
and  minimized  the  non-zero  terms  in  the  global  stiffness 
matrix. 


The  procedure  used  to  solve  the  load  carrying  weld 
geometry  is  summarized  here.  The  solution  of  other  crack 
geometries  was  similar.  The  first  step  was  to  divide  the 
structure  into  sub-structures.  A  total  of  19  substructures 
were  used.  The  stiffness  matrix  [A]  for  a  substructure  can 
be  expressed  as 


Fa  .  _  :  Aaa  ^0  ^  6a 

F0  L  A0a  A00i  .V 


(Al) 


where  subscript  a  refers  to  the  exterior  connected  nodes  of 
the  substructure  and  0  the  interior  nodes.  If  no  loads  are 
applied  to  the  interior  nodes,  this  reduces  to 


Fq  =  0  = 


A0a 


+  A 


00  U0 


(A2) 


The  displacement  of  the  interior  nodes  can  be  determined  as 


6 


0 


Afl„ 

0a  a 


(A3) 


Substituting  Eq.  A3  into  Eq.  Al  results  in  the  reduced 
fully  populated  stiffness  matrix  of  the  substructure. 
Equation  A4  relates  the  exterior  nodal  forces  and  displace¬ 


ments. 
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Fa  “  [Aaa  "  Aag  Agg  1  Aga]  6a 


(A4) 


Symmetry  of  the  stiffness  matrix  requires  that 


A  ga  '  Aag 


(A5) 


Hence,  Eq.  A4  can  be  expressed  as: 


Fa  [Aaa  "  Aag  Agg  Aag]  6a 


(A6 ) 


Only  Aaa ,  Aag ,  and  A^  need  be  generated  in  the  program. 


The  displacement  boundary  conditions  along  the 
planes  of  symmetry  were  accounted  for  by  multiplying  the 
diagonal  terms  associated  with  these  displacements  by  one 
million.  This  is  analogous  to  considering  the  displacement 
boundary  conditions  to  be  enforced  by  stiff  one  way  con¬ 
nections  or  having  a  stiffness  of 


AROLLER 


A..  [106-1] 

11 


where  A^  is  the  value  of  the  diagonal  terms  associated 
with  the  displacement  boundary  condition.  The  consequence 
of  this  modification  to  the  stiffness  matrix  is  that  the 
boundary  node  becomes  an  interior  node  connected  to  the 
stiff  roller.  Therefore  the  nodes  along  the  lines  of 
symmetry  can  be  treated  as  interior  nodes  and  included  in  the 
g  listing. 


r 
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After  the  reduced  stiffness  matrix  for  each  sub¬ 
structure  is  found  they  were  combined.  For  example,  consider 
two  substructures  i  and  j  with  reduced  stiffness  matrice  A^ 

and  A.  and  connected  at  nodes  v.  The  other  unconnected 
1 

nodes  are  a  for  substructure  i  and  1 1  for  substructure  j. 
Equilibrium  of  the  nodes  y  requires  that 

F  +  F  =  0  (A7 ) 

i  Yj 


and  compatibility  of  the  displacements  requires  that 


(A8) 


Partioning  the  reduced  substructures  stiffness  matrices 
yields  the  stiffness  matrix  for  the  combined  substructures 
as 


A 


ay 


A 


YY 


(A9) 


and 


'  Fr 


3v 

!aBy 


A 


Y0 


The  solution  of  Eq.  A7-A10  yields 


"a  _  Aaa  Aay  ^Aya 

"  -Au  FA  „ 

0  ;  0Y  Ya 


-A„  FA  a  ! 
o.y  yd 

Au  2  -  A0  FA  Q  I 
S3  0y  y0  J 


where 


YY 


+  A  ] 
YY  j 


-1 


a 


(A10) 


(All) 


F  =  (A. 


(A12) 
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The  substructures  to  the  left  of  the  crack  were  combined 
using  Eq.  All.  This  yielded  the  stiffness  matrix  for  this 
part  of  the  structure.  This  was  repeated  to  the  right  of  the 
crack.  These  stiffness  matrices  are  independent  of  crack 
length  and  consequently  are  only  determined  once  for  each 
joint  geometry  analyzed. 

The  stiffness  matrices  of  the  substructures  near 
the  crack  were  determined  for  each  crack  length.  They  were 
then  combined  with  the  stiffness  matrices  for  the  areas  to 
the  right  and  left  to  yield  the  stiffness  matrix  relating 
the  displacement  at  the  loaded  nodes  to  their  nodal  forces. 
This  12  x  12  matrix  was  inverted  and  multiplied  by  the  known 
nodal  loads  to  determine  the  nodal  point  displacements. 


The  size  of  the  matrices  to  be  inverted  by  this 

solution  scheme  are  kept  to  a  minimum.  The  inversions 

required  are  dependent  on  the  number  of  interior  nodes, 

matrix  AQO  in  Eq.  A4 ,  and  the  number  of  connected  nodes, 

PP 

matrix  F  ^  =  A  +  A  ,  in  Eq.  All.  The  time  required  for 
YVj  VYj 

the  symmetric  inversion  routine  used  was  proportional  to  the 


cube  of  the  matrix  size. 


By  dividing  the  structure  into 


small  substructures  the  time  required  for  inversion  was 


minimized.  As  the  size  of  the  matrix  A^^  decreased,  fewer 
non-zero  off  diagonal  terms  were  considered. 
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Another  advantage  of  the  method  of  substructures 
was  that  the  stiffness  matrix  of  the  segment  of  the  main 
plate  from  the  toe  of  the  weld  to  the  loaded  nodes  did 
not  have  to  be  calculated  for  each  joint  geometry.  This 
segment  geometry  is  only  a  function  of  the  main  plate 
thickness,  tp.  The  stiffness  matrix  of  this  segment  was 
independent  of  tp.  Therefore  the  stiffness  matrix  of  this 
segment  need  only  be  developed  once  for  all  joint  geometries. 
This  allows  a  considerable  reduction  in  solution  time. 

The  programs  used  for  both  the  compliance  and 
stress  analysis  made  use  of  mesh  generation  subroutines 
which  determined  the  nodal  point  locations  and  element 
coordinates.  This  minimized  the  error  associated  with  the 
generation  of  the  element  stiffness  matrix.  All  calculations 
were  performed  on  a  CDC  6400  in  single  precision. 


r 
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